iiewwY^ 


NWCroSTGWOUATt  SCHOOl 
MONTERfY  CAIIFORNW  93943 


NAVENVPREDRSCHFAC 
CONTRACTOR  REPORT 
CR  85-05 


SHIP-TO-SHIP  PROPAGATION: 
MARITIME  DATA  ANALYSES 


•  1^2)  )S5^ 


Joseph  L.  Manning,  James  A.  Dowling,  John  R.  Hummel 

OptiMetrics,  Inc. 

Ann  Arbor,  Ml  48104 


Contract  No.  N00228-82-C-8171 


MARCH  1985 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED 


Prepared  For: 

NAVAL  ENVIRONMENTAL  PREDICTION  RESEARCH  FACILITY 
MONTEREY,  CALIFORNIA  93943-5006 


K 


aUALIFIED  REQUESTORS  MAY  OBTAIN  ADDITIONAL  COPIES 
FROM  THE  DEFENSE  TECHNICAL  INFORMATION  CENTER. 

ALL  OTHERS  SHOULD  APPLY  TO  THE  NATIONAL  TECHNICAL 
INFORMATION  SERVICE. 


M  Ph 

1  ) 

A  D  ~~  w  1 5  h*  4  b  4 

1-  0 

C2) 

r’  b' 

^  CD 

C  .i 

■  •  J  )  . 

• '  u  ) 

M 

P' 

jr  .  X  lib.  ;  »■•(  X  Cb  x  CC 

*i  1 

'  w  ■• 

bn  i  D-t  o-bh D  Fr  ooaq  at on  1 

♦  c 

•  a ) 

( U ) 

I.>r4 

r  9  1 

F'j.  nax  rept. «  « 

.  1  i;j 

itanni  nq  ,  ■  j  ,  »... , 

AU 

\  1  a) } 

0 Vv 4  i r i Cl « c‘  p  He 

Hu 

S  i0y 

r‘tHfriiTi0x  n  M  r’\  M 

AD 

'•  ^  1  '■ 

hiar  i.  R'ou 

r'0 

1 C ; 

•i.  .2  /  0 

Rb 

^  i  *+  > 

LlAl;qlx6 

1  •'  T 

U  i 

^  15> 

N 2 2 d -- d  2 •■■“  (3  .i,  7  i 

rn 

■•.id) 

Nbr'KF  “~CR-“*d5“~05 

RC 

•.  2i^)  > 

L  n  c:  .L  a  ';::v  s  x  t  x  0  c;  r  '000  r"  t: 

D£ 

\  2l-I ) 

•♦'•'.Jr  1  i  t.JAL.  KADAI-^  n  u.  1  US’";  1  :  rvHfx- 

c.  *4  'y'  i  R  D  T'j  Pi  b » M  i  3  ,  M  b  b  G  R  F  7 ...  u  It ,  i 
! J  H  £  hi  1  L-  Hu.  H  3  b.  R  b  C  G  h  R  U  i  H  1'  .l 
r  b u  1.) R  i  1) b S  p  I  i''*! H b b 3  ,  LA b r. K b 
r'  R  Q  »•“  i  L  b.  3  4  P  R  iJ  F'  A  a'  A 1  i  U  N  ^  F'  L'  i~. 

I  U  ib  H  i  P  ,  !C-i  b  R  M  P-\  u-  S  L  G  G  ri  i.  *3 

DC 

(24) 

( u  > 

ID 

(25  i 

DIAL  (.1)1  i-ferenti  ai  APsor'D  tx 
F'E.b3  /'‘j4N 

IC 

( 26 ) 

(Li.) 

AB 

(27  > 

A  i“i  3.  n  a  i  v*  SIS  D  t  t  (C  s  0  c  s  0 1 
support 1 nq  meteor  Ox  oqi cai 

.1  iilt^  L  1  c'.  A.  V  ■ 


:.  r\  £1  ■:!' 


ta  has  Deen  pertormec-  , he 
T*  i  r  St:  d  at  a  se  t  vv  a  s  t:  a.  s  /  a  ri  a  i  f  u  o  r  ri  0  cl  1 1  i-  er'  E?n  c  i  a  x 
ahsorpt  1  cji ’i  i  1  cia.i'’*'  ( i-  i.  Al.  .'  v:  ’/st  0*“'  n  r  cjrn  t n:  s  ciat a  se  c. « 


V  e  Y  "  1 1  c  <a  i  v'l^  a  t  e  r  v  a  p  C'  r ' 

Lien 

•!"•  X  t  V 

pr"  o-t  1  .L  es 

were  e:-; tracts 

Cl 

ana  i..(sed  t  o  c  a  a.  c  u,  x  a  c  0 

t  r 

anc-m 

4  ttances 

1 0  r  c.  w  a  !"■(  0  D  i 

s  e  c: 

LF  laser's,  ■  nerma...  Pi 

CUMi 

X  n  q 

c:  a ...  c  IX  .L  a  1 1 

0  n  a  1  •:“>  0  w  e  0 

made 

to  cofripat^'0  the  per*Y or 

man 

C  0  L' 

■f  a  .;.  a.  ser‘ 

C'lO  1 1  r.  1  Z  0  0  i...'  r>  1 

1"'*  Ci  a 

t  a  n  d  a r  cl  a  t  ■  c  p  r  1 0  r  e  w a  c  0 

.'•■"  p  r'“ 

dt  1  .l  a  wi  t 

ct  1  a  1'  e  r 

0  D  t  .i.  (T|  ice  Ci  U.  S  .  i"i  C!  L>  1  P .....  ••*•• 

c.;  e 

4  y'  e  0 

or  1  es 

r  ;  i'”!  e  Ci  i  'f  Y'  e  *"  0  f:- 

c  e 

a s  "t 0 u. i"( 0  to  P e  iXi a .1.  i 

^  G 

r” 

e  CO n q i t X 

C'  f‘"i  ''a  u  '3  e  u  i’“i  e  '•  0  n 

T'ne 

••;■>  0  r:  0  n  a  -.  1  a.  t  a  s  e  t  c  c.'  n  c- .: 

S'l" 

C  Or 

.1.  OP' 'C  a  1 

n  at. mospi"; er  2  c. 

^ )■  ’  a  1 '  ir>  r<''  1.  t:  1  '1  '.1 0  a  \  •  ■  ..1  ■'  --.i:*  c.! 

.1.  .C  I- 

C'  i'  f  G  e 

^T'l ea ‘ii Li r  s  i ( te ri  c  s  n  -.-.■x q  0  ;i  i'"! 

;:upDOt"  T.‘  uT  :  Di'Xj  ••;  oc 

C-.  . 

;  e  1’"^  G 

r  ’"  i  ;  ir,  a  v.'  .  T’ 

0  0 emo-”' St r  a":  .i  0 

•....'"•  L..  .1  .'U .'  •!  1^! e -la  1.- i..'. r  e  ..  U'*  ■:•<.* 

iTi  X 

t  a 

■.;e-s  vven- 

conn  area  .ui  c.i 

1... !  '“'vAiH  3  L'  a.!  .1.  '•..  :•* 

'  .  V’."' 

c:.  a ..  •::: 

•.. '  .1.  •■.:•.  •...  ...  ’.  ■  •  .  .  r 

.“  '.liC'C  C  Ct"  00^10!“ 

•' 

a  s  C'  ij  i'  a  X  i "  i  e  c  .i.  i*.  -  *1  e 

'  •"• 

■  I'  .1.  L  /■■ 

•!j  '  h  G'  l  e  1' '  •  .: ' 

.  i  ■  Cl  0  M  L.  0  it'  P  r  ■  .L  3  c: 

'iC  3 

•v  0‘~  '■"!  e  b'-“  *..  2  *•!  i.  c  r‘  Oiiic- 

•“■r 

■*  r.** '.  .,•  1 

L  *  i  •j:/.avci 

•a:-  i  c  ..k  .1  a  t  e  ■.■. 

i’  ■  a  J  SiJi  I  C  ..  c  i"'  C  . 
v  a  .u  L  '  0  V.i/  n  ^  -.0  v  WO"' 


■'.J’J. 

-  "i  C.  .1 


J 


,  >ar. 


c  .•■■■ari  sii;  i.  V 

a  c iTiosp  "i 0r  .1.  c 


Cil  (  a  :;:' 


■:3  p  a.  - j  a  ! 


h-»l. 


^  xJa .  ' 


ADA;i  ;w<X  >;  j 


•  aQ  ^ 


1  V) 


^hY  CuNiHII'm  c./r  IJK  .  L  ijr-4  •.  rtUu  iu?M  i  m 

rl  I  U F  .1  C H b.  H K h  H Ij u b*'  u  1 N  i  i'~: b,  o c. i\ b h*. H L.  M  x  C H 0 r  u K PI b  K hi 


:.9;; 


i  7 


# 


CC  kZ'd) 


UNCLASSIFIED 

SECURITY  CLASSlf^lCATlON  OF  THIS  PAOt 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


lb  RESTRICTIVE  .VIARKINGS 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


2b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


3  DISTRIBUTION /availability  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

OMI-116 


5,  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

CR  85-05 


6a.  NAME  OF  PERFORMING  ORGANIZATION 


^..Opti Metrics,  Inc. 


6b  OFFICE  SYMBOL 
(If  applicable) 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Naval  Environmental  Prediction 
Research  Facility 


6c.  ADDRESS  {City,  State,  and  ZIP  Code) 

2000  Hogback  Road,  Suite  3 
Ann  Arbor,  MI  48104 


7b.  ADDRESS  (Oty,  Sfate,  and  ZIP  Code) 

Monterey,  CA  93943-5006 


8a.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 

Naval  Air  Systems  Command 


8b  OFFICE  SYMBOL 
(If  applicable) 

(AIR-330) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

N00228-82-C-8171 


8c.  ADDRESS  (C/ty,  State,  and  ZIP  Code) 

Department  of  the  Navy 
Washington,  DC  20361 


10.  SOURCE  OF  FUNDING  NUMBERS 


program 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO 

NO 

NO. 

ACCESSION  NO 

63754N 

90182-AA 

DN993603 

11  TITLE  (Include  Security  Classification) 


Ship-To-Ship  Propagation:  Maritime  Data  Analyses 


12  PERSONAL  AUTHOR(S) 


Manning,  Joseph  L.;  Dowling,  James  Ao;  Hummel,  John  R, 


13a.  TYPE  OF  REPORT 

Final 


13b.  TIME  COVERED 
FROM  TO 


14.  {Ye^r,  Month,  Day) 


15  PAQ^^(^UNT 


16.  SUPPLEMENTARY  NOTATION 


17  COSATI  CODES 

FIELD 

GROUP 

SUB-GROUP 

04 

01 

04 

02 

18.  SUBJECT  TERMS  {Continue  on  reverse  if  necessary  and  identify  by  block  number) 

DIAL  DF  laser  atmospheric  propagation 

Atmospheric  transmission 
Water  vapor  profiles 


Long  focal  length  imaging 


19  ABSTRACT  {Continue  on  reverse  if  necessary  and  identify  by  block  number) 

An  analysis  of  two  sets  of  atmospheric  optical  and  supporting  meteorological  data  has  been 
performed.  The  first  data  set  was  taken  by  an  airborne  differential  absorption  lidar 
(DIAL)  system.  From  this  data  set,  vertical  water  vapor  density  profiles  were  extracted 
and  used  to  calculate  transmi ttances  for  cw  and  pulsed  DF  lasers.  Thermal  blooming  calcu¬ 
lations  also  were  made  to  compare  the  performance  of  a  laser  optimized  using  a  standard 
atmosphere  water  vapor  profile  with  a  laser  optimized  using  DIAL-derived  profiles.  The 
difference  was  found  to  be  small  for  the  conditions  used  here. 

The  second  data  set  consisted  of  long-path  atmospheric  transmittance  and  radiosonde 
measurements  made  in  support  of  a  long  focal  length  imaging  demonstration  (LFLID). 

Measured  transmi ttances  were  compared  with  LOWTRAN  6  transmittance  calculations.  Good 
agreement  was  obtained  in  the  3-5  micrometer  range.  Comparisons  for  the  8-12  micrometer 
region  showed  calculated  transmi ttances  about  a  factor  of  10  lower  than  measured  values. 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 

®  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT  □  DTIC  USERS 

21,  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Haggerty,  J. 

22b  TELEPHONE  (Include  Area  Code) 

(408)  646-2805 

22c.  OFFICE  SYMBOL 

NEPRF  WU  6,2-21 

DD  FORM  1473, 84  MAR 


83  APR  edition  may  be  used  until  exhausted. 
All  other  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


UNCLASSIFIED 


TABLE  OF  CONTENTS 


1.  INTRODUCTION . 1 

1.1.  Background  and  General  Description 

of  this  Study . 1 

1.2.  The  NASA  DIAL  System . 3 

2.  DESCRIPTION  OF  THE  DATA  SETS . 4 

2.1.  The  NASA  DIAL  Data  Set . 4 

2.1.1.  Introduction . 

2.1.2.  Meteorological  Conditions 

2.1.3.  The  NASA  DIAL  H2O  Data... 

2.2.  The  LFID  Support  Data  Set . 9 

2.2.1.  Meteorological  Conditions  Occurring 
During  Collection  of  the  LFLID 

Support  Data . 23 

2.2.2.  The  LFLID  Support  Long-Path 

Transmittance  Data . 29 

3.  REDUCTION  AND  ANALYSIS  OF  THE  NASA  DIAL  H2O  DATA . 36 

3.1.  Average  H2O  Profile  Over  All  Flights . 36 

3.2.  Total  Column  Water  Amounts . 40 

3.3.  Fitting  the  Water  Vapor  Data . 52 

4.  INFLUENCE  OF  H^O  PROFILES  DERIVED  FROM  THE 

DIAL  DATA  ON  DF  LASER  PROPAGATION . 69 

4.1.  Model  for  Calculation  of  DF  Laser 

Transmi  ttance . 69 

4.1.1.  Transmittance  for  a  Laser 

Power  Spectrum . 69 

4.1.2.  Transmittance  Calculation  for 

the  LIDAR  Data . 7  0 


to  ►U  ►U 


4.2.  Results  of  Transmittance  Calculations . 74 

4.2.1.  CW  Laser  Transmittance  for 

Individual  Profiles . 75 

4.2.2.  Pulsed  Laser  Transmi ttances  for 

Individual  Profiles . 85 

4.2.3.  Transmi ttances  for  Average  and 

Extreme  Profiles . 92 

4.2.4.  Transmittance  Calculations  for  a 

Standard  Atmosphere . 94 

4.3.  Description  of  Thermal  Blooming  Calculat ions . . . . 94 

4.4.  Results  of  Thermal  Blooming  Calculations . 98 

5.  REDUCTION  AND  ANALYSIS  OF  THE  LFLID  SUPPORT  AND 

RADIOSONDE  TRANSMITTANCE  DATA . 101 

6.  SUMMARY  AND  CONCLUSIONS . 114 

6.1.  The  NASA  DIAL  Data  Base . 114 

6.1.1.  Summary . 114 

6.1.2.  Conclusions  and  Recommendations . 114 

6.2.  The  LFLID  Support  Data  Base . 116 

REFERENCES . 119 

ACKNOWLEDGMENTS  . 121 

DISTRIBUTION  . 123 


11 


1 

INTRODUCTION 


1.1  BACKGROUND  AND  GENERAL  DESCRIPTION  OF  THIS  STUDY 

A  project  to  measure  DF  laser  transmittance  and  high- 
resolution  transmittance  over  pathlengths  of  several  kilo¬ 
meters  between  two  ships  was  planned  by  the  Navy  for  1982. 
These  measurements,  along  with  supporting  data  such  as  me¬ 
teorological  measurements,  would  have  generated  a  data  base 
for  analysis  and  prediction  of  problems  in  DF  laser  propaga¬ 
tion,  especially  at  sea.  Among  the  problems  which  might 
have  been  addressed  are  accuracy  of  line  parameters  (both 
for  attenuation  calculations  and  for  thermal  blooming  pre¬ 
diction);  aerosol  extinction;  validation  of  aerosol  models 
and  applicability  of  existing  models  for  open  ocean;  HDO/H2O 
ratio;  and  H2O  continuum  [1]. 

The  ship-to-ship  measurement  project  did  not  take 
place.  In  this  report  we  present  a  study  based  on  the 
analysis  of  two  backup  data  sets.  The  first  one  considered 
consists  of  water  vapor  profiles  over  open  ocean  measured 
with  an  airborne  DIAL  system.  The  second  data  set  consists 
of  slant-path  transmittance  and  radiosonde  data  measured  in 
conjunction  with  a  Long  Focal  Length  Imaging  Demonstration 
(LFLID).  These  data,  supporting  tests  conducted  by  the 
Naval  Air  Development  Center  (NADC),  were  measured  using  a 
66  km  slant-path  between  the  islands  of  Maui  and  Lanai  in 
Hawaii.  With  these  data  sets  it  was,  of  course,  not  possible 
to  address  most  of  the  problems  for  DF  propagation  mentioned 
in  the  previous  paragraph,  although  this  study  examines  the 
effects  of  measured  water  vapor  profiles  on  the  propagation 
of  high-energy  DF  laser  beams,  using  the  LIDAR  data  set. 
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The  DIAL  data  were  examined  for  realism  and  consist¬ 
ency.  Transmi ttances  were  calculated  for  pulsed  and  cw  DF 
lasers  for  vertical  and  slant  ranges  for  each  profile.  The 
grid  over  which  the  data  were  taken  included  the  boundary  of 
the  gulf  stream,  so  the  influences  of  sea  surface  tempera¬ 
ture  in  general  and  the  gulf  stream  boundary  in  particular 
were  examined.  Average  and  extreme  water  vapor  profiles 
based  on  the  data  were  generated,  and  DF  laser  transmit- 
tances  for  these  profiles  were  compared  to  transmittances 
calculated  using  standard  atmosphere  water  vapor  profiles. 

Thermal  blooming  calculations  were  done  to  determine 
whether  there  would  be  any  real  advantage  to  having  a  ship¬ 
board  device  to  monitor  the  water  vapor  profile.  The  ques¬ 
tion  we  set  out  to  answer  was  this:  if  a  laser  system  is 
optimized  for  atmospheric  conditions  based  on  surface  mea¬ 
surements  and  standard  atmosphere  vertical  profiles,  how 
much  less  effective  will  it  be  than  if  actual  vertical  mea¬ 
surements  were  used? 

The  long-path  transmittance  data  from  the  Hawaii  tests 
were  obtained  from  the  AVCO  Corporation  which  was  primarily 
responsible  for  carrying  the  transmittance  measurements  as 
part  of  their  support  of  the  LFLID  tests.  OptiMetrics,  Inc. 
(OMI)  participated  in  the  tests  by  providing  and  operating  a 
searchl ight-blackbody  source  system  on  the  island  of  Lanai 
during  the  tests.  The  transmittance  data  were  collected 
using  a  66  km  path  between  Lanai  and  Mount  Haleakala,  Maui. 
The  data  were  compared  with  predicted  transmittances,  based 
water  vapor  density  profiles  obtained  from 
The  transmittance  predictions  were  made  using 
model  [2] . 


on  vertical 
rad iosondes . 
the  LOWTRAN 
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1.2  THE  NASA  DIAL  SYSTEM 


The  water  vapor  profiles  were  measured  using  an  air¬ 
borne  DIAL  (differential  absorption  lidar)  system  developed 
by  NASA  [3] .  This  system  uses  two  high-conversion-efficiency 
tunable  dye  lasers  pumped  by  two  frequency-doubled  Nd : YAG 
lasers.  One  laser  operates  on  an  H2O  absorption-line  center 
wavelength  while  the  other  is  tuned  off  of  the  absorption 
line.  The  concentration  of  H2O  at  each  range  is  determined 
frcTO  the  strengths  of  the  returns  for  the  two  wavelengths  at 
that  range.  The  return  at  the  off-line  wavelength  can  also 
be  used  for  aerosol  measurement.  The  entire  system  is  in¬ 
stalled  in  an  Electra  aircraft. 

Data  were  taken  on  14,  18,  and  24  September  1982.  Only 
the  data  from  24  September  were  supplied  for  analysis. 
Aerosol  information  was  not  derived  fron  the  data,  and  the 
emphasis  in  our  analysis  was  on  water  vapor. 
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2 

DESCRIPTION  OF  THE  DATA  SETS 

2.1  THE  NASA  DIAL  DATA  SET 

2.1.1  INTRODUCTION 

The  DIAL  data  to  be  discussed  here  were  taken’  on  24 
September  1982  off  the  coast  of  the  United  States  east  of 
Virginia  and  North  Carolina.  The  flights  of  the  airborne 
instrument  took  place  between  1126  and  1512  local  time.  A 
total  of  eleven  flights  took  place.  Figure  1  displays  the 
area  of  the  flights  and  shows  the  flight  paths  of  the  vari¬ 
ous  flights.  The  flights  were  generally  north-south  and 
took  place  over  the  Gulf  Stream.  Table  1  summarizes  the 
flights.  The  complete  data  base  to  be  analyzed  consists  of 
meteorological  data  that  were  collected  separately  and  the 
NASA  DIAL  data. 

2.1.2  METEOROLOGICAL  CONDITIONS 

The  weather  during  the  flights  was  generally  good.  The 
weather  situation  was  one  of  clearing  and  generally  warmer 
weather  after  a  week  of  rain.  Figure  2  displays  the  surface 
weather  chart  for  1800  GMT  (1300  local  time)  for  the  United 
States.  The  test  area  was  under  the  influence  of  a  high 
pressure  system  just  to  the  north  of  the  test  site.  Table  2 
lists  the  surface  observations  from  Wallops  Island,  Virginia 
and  ships  that  were  in  the  vicinity  of  the  test  area. 
Wallops  Island  reported  scattered  clouds  at  the  times  of 
their  0000  GMT  and  1200  GMT  radiosonde  launches.  The  ships 
reported  low  clouds  with  total  amounts  over  50%.  The  values 
reported  are  consistent  with  mean  values  given  in  the  U.S. 
Navy  Marine  Climatic  Atlas  of  the  World  [4]. 
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FIGURE  1.  FLIGHT  PATHS  OF  THE  NASA-LANGLEY  AIRCRAFT 

AND  POSITIONS  OF  EXPENDABLE  BATHYTHERMOGRAPHS 
(XBT)  USED  TO  OBTAIN  OCEAN  TEMPERATURES  AND 
SHIPS  THAT  REPORTED  METEOROLOGICAL  OBSERVA¬ 
TIONS  . 
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TABLE  1.  SUMMARY  OF  FLIGHT  PATHS 


FLIGHT 

FLIGHT  TIME 
( LOCAL) 

LATITUDE  COVERAGE 
(N) 

1 

11:26:18  - 

11:47:19 

36°43'  -  35 

°58'  -  36°33’ 

2 

12:02:16  - 

12:15:36 

36°26 ’ 

-  35°37’ 

3 

12:22:12  - 

12:35:10 

35°38' 

-  36°29’ 

4 

12:59:24  - 

13 : 18 :01 

36°47 ' 

-  35°37' 

5 

13:27:48  - 

13:42:28 

35°30’ 

-  36°26’ 

6 

13:48:10  - 

14:11:55 

36°27  ' 

-  35°38’ 

7 

14:11:55  - 

14:23:38 

35°43' 

-  36°29’ 

8 

14:29:49  - 

14 : 38 :09 

36°29  ' 

-  35°59' 

9 

14:41:15  - 

14 : 50:08 

35°56' 

-  36°29' 

10 

14:51:19  - 

14:59: 59 

36°34' 

-  37°  8’ 

11 

15:  4:41- 

15:12:43 

37°26’ 

-  37°48' 
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TABLE  2 


SURFACE  METEOROLOGICAL  OBSERVATIONS 
FOR  24  SEPTEMBER  1982. 


Wallops  Island , 

Virginia  (37°51' 

'N,  75°29'W) 

TIME 

(GMT) 

SURFACE 

TEMPERATURE 

(C) 

RELATIVE 

HUMIDITY 

(%) 

WIND 

SPEED 
(KNOTS ) 

WIND 

DIRECTION 

(DEGREES) 

1200 

10 . 6 

.86 

1 

350 

0000 

18.6 

63 

5 

155 

Ship  Observations 


TIME 

(GMT) 

TEMPERATURE 

AIR  OCEAN 

(C) 

RELATIVE 

HUMIDITY 

(%) 

WIND 

SPEED 

(KNOTS) 

WIND 

DIRECTION 
( DEGREES) 

0600j 

21.7 

24.4 

56 

9 

20 

0600^ 

19.9 

26.7 

52 

5 

320 

1200^ 

18.8 

12.7 

78 

8 

350 

1800^ 

21 . 1 

27 . 2 

62 

8 

40 

1800^ 

22.0 

— 

67 

6 

40 

1800° 

21 . 2 

21.0 

60 

9 

90 

Ship  Locations: 


1 

2 

3 

4 

5 

6 


36 .2°N, 
36.9°N, 
36.8°N, 
35.9°N, 
36.5°N, 
37 .5°N, 


74 . 5°W 
74 . 6°W 
74 . 7°W 
74 . 9°W 
74 . 2°W 
74 . 6°W 
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Two  instrumented  buoys  were  operating  in  the  vicinity 
of  the  test  area  and  obtained  meteorological  data  on  the 
days  of  the  flights.  Those  data  are  summarized  in  Table  3. 
Finally,  data  were  obtained  by  two  ship-dropped  expendable 
bathythermographs  (XBT) .  The  XBT ' s  measure  water  tempera¬ 
ture  as  a  function  of  depth.  Those  data  are  given  in  Table 
4. 

2.1.3  THE  NASA  DIAL  H2O  DATA 

The  data  obtained  by  the  instrument  consisted  of  water 
vapor  volume  mixing  ratios  as  a  function  of  altitude  in  a 
vertical  column.  The  data  were  given  at  15  m  intervals. 
The  length  of  the  vertical  column  varied  from  one  location 
to  another  but  covered  the  range  131  to  1361  m.  The  instru¬ 
ment  took  data  about  every  20  seconds.  Figures  3  through  13 
give  the  average  water  vapor  profile  for  each  of  the 
flights.  There  is  a  fair  amount  of  structure  in  the  profile 
but  they  generally  show  decreasing  values  with  altitude. 
The  dramatic  changes  in  the  profiles  at  the  top'  and  bottom 
of  the  profiles  are  believed  to  be  indications  of  the  limits 
of  sensitivity  of  the  instrument  and  not  physical.  The 
magnitudes  are  in  general  agreement  with  climatological 
averages  for  that  latitude  [5]. 

2.2  THE  LFID  SUPPORT  DATA  SET 

The  meteorological  conditions  occurring  during  the  col¬ 
lection  of  the  long-path  transmittance  measurements  were 
fairly  constant.  The  most  useful  long-path  measurements 
were  made  during  the  early  morning  hours  on  4  days  in  July 
1984.  Although  long-path  transmittance  and  radiosonde 
measurements  were  collected  on  several  additional  days,  on 
these  occasions  the  long-path  transmittance  measurements 
were  beset  by  and  at  least  partially  contaminated  by  the 


9 


TABLE  3.  METEOROLOGICAL  DATA  OBTAINED  FROM  INSTRUMENTED 
BUOYS  NEAR  THE  TEST  AREA. 


BUOY:  41001 

LATITUDE 

3^  4  1*  59** 

N  LONGITUDE 

if  18'  oo" 

W 

WIND 

SEA  SURFACE 

TIME 

BOTTOM  DEPTH 

AIR  TEMP 

PRESSURE 

WINDSPEED 

DIRECTION 

TEMP 

(GMT) 

(m) 

(C) 

( mb ) 

( m/s ) 

( degrees ) 

(C) 

0 

0 

4353 

22 . 7 

1018  9 

5.2 

308 .0 

25.02 

1 

0 

4353 

22.7 

1019.7 

3.8 

315.5 

24 .98 

2 

0 

4353 

22 .6 

1019.8 

4  .  1 

322  1 

24  93 

3 

0 

4353 

22.5 

1020. 1 

3.7 

329.8 

24 .92 

4 

0 

4353 

22.2 

1020. 1 

2.3 

305  7 

24 .88 

5 

0 

4353 

22  .  2 

1019.9 

2 .8 

306 . 5 

24 . 84 

6 

0 

4353 

22.  1 

1020.0 

3.2 

307  .0 

24 .88 

7 

0 

4353 

22 .0 

1019.9 

2.6 

329 . 3 

24 . 87 

8 

0 

4353 

21.9 

1019  6 

3.6 

324.3 

24.87 

9 

0 

4353 

22 .0 

1019.3 

3.8 

314.7 

24  86 

10 

0 

4353 

22 .0 

1019.3 

4 .7 

315.6 

24  86 

1  1 

0 

4353 

22.0 

1019.6 

5.7 

348.2 

24  85 

12 

0 

4353 

22 .0 

1019.9 

6 . 2 

27 .0 

24  84 

13 

0 

4353 

22.0 

1020.0 

6.8 

33.0 

24 .86 

14 

0 

4353 

22.  1 

1020.3 

6 . 3 

49.0 

24  88 

15 

0 

4353 

22 . 2 

1019.9 

7  .  2 

52 . 2 

24 .89 

16 

0 

4353 

22 . 5 

1019.5 

5.2 

47.2 

24.92 

17 

0 

4353 

22  6 

1019.3 

5.6 

50. 3 

24 . 94 

18 

0 

4353 

22  7 

1018.7 

5.3 

49.2 

24 .96 

19 

0 

4353 

22 . 8 

1017 .9 

5.0 

42 . 2 

24 .99 

20 

0 

4353 

22 .8 

1016.9. 

6.4 

37 . 3 

25.00 

21 

0 

4353 

22 . 9 

1017.3 

7  .0 

26  7 

24 .99 

22 

0 

4353 

23  0 

1017 . 1 

8.5 

29.0 

24 .98 

23 

0 

4353 

23  .  1 

1016.7 

9.0  . 

32 . 2 

24  96 

BUOY:  41002 

LATITUDE 

320  18*  00"N  LONGITUDE 

7^  18*  OCT' 

W 

WIND 

SEA  surface 

TIME 

BOTTOM  DEPTH 

AIR  TEMP 

PRESSURE 

WINDSPEED 

DIRECTION 

TEMP 

(GMT  ) 

(m) 

(C) 

( mb) 

(  m/s  ) 

( degrees ) 

<C) 

0 

0 

3749 

24  6 

1018 .5 

6 . 2 

17  2 

28.00 

1 

0 

3749 

24 . 7 

1019  0 

5.7 

46.  1 

27 .97 

2 

0 

3749 

24  7 

1019.2 

6 . 2 

38 . 3 

27 . 95 

3 

0 

3749 

24 . 6 

1019.4 

6.2 

40 . 4 

27 . 95 

4 

0 

3749 

24 .5 

1019.5 

5.9 

47 . 5 

27 . 95 

5 

0 

3749 

24 .5 

1019.4 

5.6 

37 . 6 

27 . 95 

6 

0 

3749 

24 .4 

1019.0 

5.4 

40.0 

27  95 

7 

0 

3749 

24 . 5 

1018  3 

5.9 

44 .0 

27 . 92 

8 

0 

3749 

24 . 4 

1018  0 

6 . 1 

37.5 

27 .92 

9 

0 

3749 

24.5 

1017.7 

6 . 8 

41.6 

27 .92 

10 

0 

3749 

24 .5 

1018.2 

6.6 

40. 2 

27 .90 

1  1 

0 

3749 

24 .4 

1017.7 

7 . 4 

36.4 

27 .90 

12 

0 

3749 

24  4 

1018.2 

6.8 

36 . 3 

27.90 

13 

0 

3749 

24.5 

1018.3 

6.9 

43.7 

27.90 

14 

0 

3749 

24 . 6 

1018 . 1 

7 .0 

43 . 2 

27 .88 

15 

0 

3749 

24.7 

1018.2 

6.8 

48 . 7 

27 . 88 

16 

0 

3749 

24 . 7 

1018.0 

6.5 

51.7 

27  88 

17 

0 

3749 

24.7 

1017 .9 

6.3 

39 . 5 

27 . 88 

18 

0 

3749 

24.7 

1017 . 1 

6.4 

47 . 4 

27 . 88 

19 

0 

3749 

24.7 

1015.8 

6.  1 

45 . 4 

27 .88 

20 

0 

3749 

24 . 8 

1015  8  ^ 

5.8 

50.4 

27.88 

21 

0 

3749 

25.0 

1015  5 

5.5 

43.  1 

27  88 

22 

0 

3749 

25.  1 

1015 . 1 

5.4 

60.1 

27  90 

23 

0 

3749 

25.2 

1014.6 

5.5 

57 .5 

27 .90 
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TABLE  4 


DATA  OBTAINED  FROM  EXPENDABLE  BATHYTHERMOGRAPHS 


Time 

Position 

Temperature 

Ocean 

(GMT) 

Latitude 

Longitude 

(C) 

Depth  (m) 

1:20  . 

37°26 'N 

75°32'W 

21.63 

0 

2:48 

37°28  'N 

75  =  22  'w 

21 . 27 

0 
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Altitude  (m) 


Flight  1 


Volume  Mixing  Rotio  (pp  10000) 


'FIGURE  3.  AVERAGE  WATER  VAPOR  PROFILE  FOR 
FLIGHT  1. 
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Altitude  (m) 


Flight  2 

— h- 


80  100  120  140  160  '  180  200 

Volume  Mixing  Ratio  (pp  10000) 


FIGURE  4.  AVERAGE  WATER  VAPOR  PROFILE  FOR 
FLIGHT  2. 
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Altitude  (m) 


Flight  3 


+ 


H - ^ - 1 - ^ - h 

60  80  100  120  140 

Volume  Mixing  Ratio  (pp  10000) 


FIGURE  5  AVERAGE  WATER  VAPOR  PROFILE  FOR 
FLIGHT  3. 
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Altitude  (m) 


Flight  4 


FIGURE  6.  average  WATER  VAPOR  PROFILE  FOR 
FLIGHT  4. 
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Altitude  (m) 


Flight  5 


100  120  140  160 

Volume  Mixing  Ratio  (pp  10000) 


FIGURE  7.  AVERAGE  WATER  VAPOR  PROFILE  FOR 
FLIGHT  5. 
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Altitude  (m) 


Flight  6 


FIGURE  8.  AVERAGE  WATER  VAPOR  PROFILE  FOR 
FLIGHT  6. 
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Altitude  (m) 


Flight  7 


FIGURE  9 .  AVERAGE  WATER  VAPOR  PROFILE  FOR 
FLIGHT  7. 
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Altitude  (m) 


Flight  8 


FIGURE  1C .  AVERAGE  WATER  VAPOR  PROFILE  FOR 
FLIGHT  8. 
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Attitude  (m) 


Flight  9 


H ^ - 1 - ^ ^ ^ - ^- 

40  60  80  100  120  140  160 

Volume  Mixing  Ratio  (pp  10000) 


FIGURE  11.  AVERAGE  WATER  VAPOR  PROFILE  FOR 
FLIGHT  9. 
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Altitude  (m) 


Flight  10 


FIGURE  12.  AVERAGE  WATER  VAPOR  PROFILE  FOR 
FLIGHT  10. 
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Altitude  (m) 


Flight  11 


FIGURE  13.  AVERAGE  WATER  VAPOR  PROFILE  FOR 
FLIGHT  11. 
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presence  of  some  degree  of  cloud  cover  along  the  measurement 
1 ine-of-s ight .  The  four  days  during  which  the  most  reliable 
transmittance  data  were  collected  were  July  6,  July  12,  July 
23,  and  July  25,  1984.  Table  5  lists  the  dates  and  times 

for  which  the  transmittance  and  corresponding  radiosonde 
data  discussed  in  this  report  were  collected.  Also  listed 
in  Table  5  are  the  designations  of  the  individual  measured 
spectra.  (See  the  note  accompanying  the  Table  for  an  explan¬ 
ation  of  these  designations  which  are  used  to  identify  the 
measured  spectra  in  the  plots  of  the  data  presented  in 
Section  2.2.2) . 

2.2.1,  METEOROLOGICAL  CONDITIONS  OCCURRING  DURING  COL¬ 
LECTION  OF  THE  LFLID  SUPPORT  DATA 

Meteorological  conditions  were  monitored  during  the 
LFLID  support  tests  by  radiosondes  launched  from  a  sea-level 
location  in  Kihei,  Maui.  This  site  is  located  approximately 
21  km  from  the  Air  Force  Maui  Optical  Station  (AMOS)  site  at 
the  top  of  Mt.  Haleakala  on  Maui.  The  receiving  telescopes 
used  in  collecting  the  long-path  transmittance  data  were 
located  at  the  AMOS  observatory.  Figure  14  shows  the  loca¬ 
tion  of  the  various  sites  used  in  the  transmittance  and 
radiosonde  measurements.  The  transmittance-measurement 

receiver  telescopes  were  located  at  an  altitude  of  3.05  km 
at  the  AMOS  observatory  indicated  by  the  symbol  R  in  the 
Figure.  The  sources  used  in  the  measurement  were  located 
approximately  40  m  above  sea  level  near  Maneli  Bay  on  the 
island  of  Lanai  (indicated  by  the  symbol  S  in  the  Figure). 
The  slant-path  range  between  the  source  and  receiver  was 
66.2  km . 

A  preliminary  set  of  measurements  was  made  in  January 
1984  with  the  source  located  at  the  position  S'  in  the 
foothills  of  the  West  Maui  mountains,  using  a  35  km  path. 
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TABLE  5.  TRANSMITTANCE  AND  RADIOSONDE  MEASUREMENT  TIMES 


TRANSMITTANCE 

RADIOSONDE 

TRANSMITTANCE 

MEASUREMENT 

MEASUREMENT 

SPECTRUM 

DATE 

DATE 

DESIGNATION 

GMT  HST* 

GMT  HST 

7-6-84 

7-6-84 

1012  0012 

1044  0044 

S06R1** 

1052  0052 

1044  0044 

L06R2 

1242  0242 

1044  0044 

S06R5 

7-12-84 

7-12-84 

1208  0208 

1311  0311 

S12R1 

7-23-84 

7-23-84 

1209  0209 

1130  0130 

S23R1 

1246  0246 

1130  0130 

L23R2 

1324  0324 

1130  0130 

S23R3 

1402  0402 

1130  0130 

L23R4 

7-25-84 

7-24/25-84 

1051  0051 

0836  2236 

(7/24) 

S25R1 

1226  0226 

0836  2236 

(7/24) 

L25R3 

1303  0303 

0836  2236 

(7/24) 

S25R4 

1343  0343 

0836  2236 

(7/24) 

L25R5 

1419  0419 

0836  2236 

(7/24) 

S25R6 

=  HST  (Hawaiian  Standard  Time) 

=  A  code  identifying  a  particular  measured  transmittance 
spectrum.  The  initial  S[L]  indicates  the  short  wave¬ 
length  (3-5  pm)  or  [long]  wavelength  (8-12  ym)  spectral 
band.  The  next  two  digits  indicate  the  date  in  July 
1984  when  the  spectrum  was  measured  and  the  last  two 
digits  indicate  the  sequential  position  of  the  particu¬ 
lar  data  file  recorded  on  that  date;  e.g.  S23R3  indi¬ 
cates  a  short  wavelength  spectrum  recorded  on  July  23 
in  the  third  sequential  file  for  that  date. 
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FIGURE  14.  SITES  USED  IN  LFLID  SUPPORT  MEASUREMENTS. 


The  data  from  the  earlier  tests  are  not  considered  here, 
however.  The  site  in  Kihei  used  for  the  radiosonde  launches 
is  shown  by  the  symbol  M  in  Figure  14. 

Figures  15  and  16  show  the  measured  temperature  and 
derived  water  vapor  density  profiles,  respectively,  for  the 
radiosonde  data  which  accompanied  the  selected  long-path 
transmittance  spectra  discussed  in  this  report. 

The  temperature  data  shown  in  Figure  15  show  the  pres¬ 
ence  of  temperature  inversions  occurring  at  altitudes  be¬ 
tween  1.5  and  2.1  km  in  three  of  the  four  measurement 

periods.  Only  the  data  for  6  July  show  a  consistent  lapse 
of  temperature  with  increasing  altitude. 

The  four  water  vapor  density  profiles  corresponding  to 
the  transmittance  measurements  are  shown  in  Figure  16.  They 
were  derived  from  radiosonde-measured  temperature  and  rela¬ 
tive  humidity  profiles.  They  show  a  rather  consistent  de¬ 
crease  in  water  vapor  density  ( PH2O)  with  a  relatively  small 
variation  at  each  altitude  (~  ±  1.3  g/m^  from  the  average 

value  at  a  given  altitude)  for  altitudes  below  2.1  km. 

Above  this  level  a  larger  spread  in  the  values  is  observed 

ranging  from  0.2  g/m^  to  5.8  g/m^  at  an  altitude  of  3  km. 
(The  AMOS  observatory  is  located  at  an  altitude  of  3.05 
km .  ) 

The  source  (S  in  Figure  14)  was  located  sufficiently 
far  from  the  ocean  shore  that  locally  high  aerosol  scat¬ 
tering,  related  to  heavy  surf  conditions,  would  not  contri¬ 
bute  to  the  measured  transmi ttances .  Sugar  cane  fires, 
generating  a  large  amount  of  smoke,  are  present  during  the 
sugar  cane  harvest  and  normally  occur  during  daytime  and 

early  evening.  The  long-path  transmittance  measurements  were 
carried  out  during  the  hours  between  midnight  and  dawn,  well 
after  smoke  from  the  burning  cane  fields  had  been  ventilated 
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ALTITUDE  (KM) 


FIGURE  15.  VERTICAL  TEMPERATURE  PROFILES  MEASURED  WITH 
A  RADIOSONDE. 
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ALTITUDE  (KM) 


FIGURE  16.  VERTICAL  WATER  VAPOR  DENSITY  PROFILES 
DERIVED  FROM  RADIOSONDE  MEASUREMENTS . 
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away  from  the  path  used  for  the  transmittance  measurements. 
Trade  winds,  a  dominant  feature  of  the  normal  weather  pat¬ 
terns  of  the  Hawaiian  islands,  typically  persist  from  a 
northeasterly  direction  at  speeds  varying  between  4  and  10 
m/s,  and  would  clear  the  transmittance-measurement  path  of 
island-generated  aerosols  during  the  few  hours  after  the 
cane  fields  had  been  burned. 

2.2.2  THE  LFLID  SUPPORT  LONG-PATH  TRANSMITTANCE  DATA 

2. 2.2.1  Measurement  Configurations  and  Instrumentation 

The  instrumentation  used  in  the  LFLID  support  transmit¬ 
tance  measurements  included  a  1200  K  blackbody  source,  used 
alternately  with  two  optical  projection  systems  for  two  sep¬ 
arate  measurement  configurations.  Two  large-aperture  re¬ 
ceiver  telescopes  were  also  used  in  the  measurements.  One 
measurement  configuration  served  as  the  basis  for  providing 
an  absolute  transmittance  calibration  for  the  long-path 
measurements.  In  this  configuration  the  radiant  output  from 
the  blackbody  source  was  projected  through  the  optical 
system  of  an  eight-inch  aperture,  f/6  Newtonian  telescope. 

For  the  absolute-transmittance-measurement  configura¬ 
tion  the  AMOS  Advanced  Multicolor  Tracker  or  AMOS  (AMTA) 
sensor  package  [6] ,  mounted  on  the  AMOS  1.2  m  aperture,  B29 
telescope,  served  as  the  receiver  system.  The  AMTA  package 
consists  of  a  25  element  infrared  sensor  array,  a  dual  chan¬ 
nel  contrast-mode  photometer  and  a  boresight  TV  camera. 

A  second  measurement  configuration  was  used  to  obtain 
moderate-resolution  (nominally  two-percent  of  the  base  wave¬ 
length)  spectral  data.  In  this  configuration  the  blackbody 
source  was  used  with  a  60-inch  aperture  optical  projection 
system  built  into  a  modified  World  War  II  anti-aircraft 
searchl ight . 
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Figure  17  shows  a  construction  diagram  for  the  modified 
60"  searchlight  which  is  used  both  as  a  source  and  receiver 
optical  system  in  long-path  atmospheric  transmittance  mea¬ 
surements.  The  major  modification  consists  in  the  use  of  a 
two-sided  secondary  mirror  to  provide  a  Cassegrain ian  opti¬ 
cal  system  when  used  with  the  original  searchlight  reflector 
(a  60"  diameter,  f/0.4,  electro-formed  parabolic  mirror). 
The  incoming  radiation  blocl<ed  by  the  15"  diameter  secondary 
mirror  is  collected  at  a  second  focal  point  formed  by  a 
parabolic  surface  incorporated  into  the  back  of  the 
secondary  mirror. 

Figure  18  contains  optical  diagrams  showing  the  geome¬ 
try  and  calculated  performance  of  the  modified  searchlight 
optical  system  when  used  with  a  one-inch  diameter  blackbody 
source  as  was  the  case  in  the  LFLID  support  measurements. 

The  searchlight  transceiver  optical  system  was  provided 
for  the  LFLID  support  tests  by  the  U.S.  Army  Atmospheric 
Sciences  Laboratory  (ASL)  at  White  Sands  Missile  Range,  New 
Mexico.  It  was  operated  during  the  tests  by  OMI  under  con¬ 
tract  to  ASL. 

The  design  and  construction  of  the  modified  optical 
system  was  originated  by  Richard  Horton  of  the  AVCO  Corpo¬ 
ration,  while  formerly  with  the  Naval  Research  Laboratory 
(NRL)  Washington,  D.C.  This  system  was  provided  to  ASL  by 
NRL  for  use  in  long-path  atmospheric  transmittance  measure¬ 
ments  . 

The  receiver  optical  system  used  when  the  60-inch 
optical  projection  system  was  in  use  was  the  AMOS  1.6  m 
aperture  telescope,  coupled  to  the  AMOS  Long  Wavelength 
Infrared  (LWIR)  Spectral  Radiometer  [6]  .  The  LWIR  sensor 
package  incorporates  a  two-element  Si:As  detector  cooled  to 
11  K,  providing  a  spectral  coverage  from  2  to  24  wm.  A 
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FIGURE  17.  60"  SEARCHLIGHT  SOURCE  CONSTRUCTION  DIAGRAM. 
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FIGURE  18.  MODIFIED  60"  SEARCHLIGHT  OPTICAL  SYSTEM  PERFORMANCE. 


spectral  resolution  of  2%  of  the  base  wavelength  in  two 
bands  (3.2  to  5.5  wm  and  7.8  to  14.0  pm)  is  provided  by  the 
use  of  two  circular  variable  filters  (CVF).  The  instan¬ 
taneous  f ield-of-view  of  the  LWIR  spectral  radiometer  is 

10x20  arcseconds. 

The  transmittance  signal  from  the  cooled  preamplifier 
in  the  LWIR  pac)cage  was  processed  with  a  Princeton  Applied 
Research  (PAR)  model  HR-8  loc)c-in  amplifier  (LIA).  The 

lock-in  amplifier  was  supplied  with  a  reference  signal  ob¬ 

tained  from  a  chopper  mounted  at  the  source.  The  chopper- 
reference  was  transmitted  to  the  AMOS  observatory  using  a 
radio  link.  A  155  MHz  single-side-band  (SSB)  FM  transmitter 
and  receiver  were  operated  between  the  two  ends  of  the 
transmittance  measurement  path  for  this  purpose.  Voice 
communications  were  available  when  the  reference  signal  was 
not  required. 

The  analog  output  signal  from  the  LIA  was  sent  to  an 

A  LI  12  bit  A/D  digitizer,  and  sampled  on  program  command 
using  a  Commodore  64  computer.  The  analog  output  from  the 
CVF  passband  wavelength  encoder  was  digitized  using  another 
channel  of  the  digitizer  for  simultaneous  sampling  by  the 
data  sampling  program. 

The  CVF  controller  was  operated  in  a  step  and  hold 
mode,  with  101  equal-wavelength  steps  covering  each  of  the 
two  wavelength  bands.  Optimally,  the  transmittance  could  be 
sampled  at  1%  wavelength  intervals,  requiring  67  samplings 
and  a  varying  wavelength  interval.  It  was  more  straight¬ 
forward  to  sample  the  data  using  a  larger  number  of  equal 
wavelength  intervals.  The  data  sampling  algorithm  allows  a 
settling  time,  approximately  equal  to  the  lock- in-ampl i f ier 
time  constant,  between  samples  to  allow  low  level  micro¬ 
phonics  produced  by  the  CVF  stepper  motor  to  decay.  The 
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measured  spectrum  can  be  observed  on  a  slowly-sweeping 
signal  analyzer. 

The  algorithm  used  for  data  collection  was  optimized  to 
allow  data  collection  under  differing  noise  conditions  and 
for  different  correlation  times  for  the  fluctuations  in  the 
transmittance  signal. 

2. 2.2. 2  Calibration  of  the  Long-Path  Transmittance  Data 

Two  types  of  calibration  measurements  were  performed  to 
provide  an  absolute-transmittance  calibration  of  the  lonq- 
path  spectra  obtained  with  the  LWIR  Spectral  Radiometer. 
The  radiometric  response  of  the  AMTA  sensor  package  was 
calibrated  using  a  "Jones  Calibration"  at  the  aperture  of 
the  1.2  m  telescope.  In  this  measurement  the  AMTA  detector 
element  is  overfilled  by  the  blackbody  radiant  output 
collected  by  the  1.2  m  telescope  and  imaged  at  the  AMTA 
detector.  No  collimating  optics  are  used  with  the  source  in 
this  configuration. 

A  second  radiometric  calibration  was  performed  using 
the  blackbody  source  together  with  the  eight-inch  Newtonian 
telescope  collimating  optics.  In  this  case  two  "short-path" 
measurement  geometries  were  used.  One  path  was  approximately 
100  m  in  length.  In  this  case  the  source  and  source  optics 
were  located  near  the  AMOS  observatory  on  the  top  of  Mt. 
Haleakala.  The  second  location  used  for  the  source  and  its 
collimating  optics  was  about  4  km  distant  from  the  AMOS 
observatory  at  a  site  near  the  observatory  access  road  which 
could  be  viewed  by  the  two  AMOS  telescopes. 

The  two  latter  measurements  provided  a  calibration  of 
the  overall  measurement  system  efficiency  including  the 
source  and  receiver  optical  systems. 
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The  long-path  spectra  obtained  with  the  60"  search- 
1 ight-blackbody  source  and  1.6  m  telescope  -  LWIR  sensor 
receiver  were  normalized  for  absolute  transmittance  by  means 
of  the  independent  measurements  made  using  the  same  IR 
source  projected  through  the  8"  Newtonian  telescope  and 
received  by  the  1.2  m  telescope  -  AMTA  sensor  combination. 
Convolution  integrals  of  the  AMTA  filter  functions  with  the 
measured  LWIR  relative  transmission  spectra  were  performed 
for  each  measured  spectrum  and  the  results  were  equated  with 
the  measured  transmittance  for  that  particular  AMTA  filter 
band.  The  resulting  multiplicative  scale  factors  thus  ob¬ 
tained  were  used  to  provide  the  absolute  normalization  of 
the  LWIR  Spectral  Radiometer  spectra.  The  estimated  uncer¬ 
tainties  in  the  absolute  calibration  of  the  resulting 
spectra  are  approximately  ±  30%. 


Prior  to  the  absolute  transmittance  normalization 
described  above,  instrumental,  wavelength-dependent  effects 
were  removed  from  the  LWIR  Spectral  Radiometer  data  by 
ratioing  each  long-path  spectrum  to  a  "zero-path"  spectrum. 
The  latter  were  obtained  using  the  blackbody  source  and  60" 
projection  optics  with  the  1.6  m  telescope  and  LWIR  Spectral 
Radiometer  operating  over  the  100  m  path  near  the  AMOS 
observatory . 
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3 

REDUCTION  AND  ANALYSIS  OF  THE  NASA  DIAL  H2O  DATA 

3.1  AVERAGE  H2O  PROFILE  OVER  ALL  FLIGHTS 

An  average  water  vapor  profile  was  obtained  from  the 
DIAL  data  by  averaging  all  of  the  data  from  flights  1 
through  10  for  each  altitude.  The  resulting  profile  is 
shown  in  Figure  19.  There  were  only  a  few  soundings  with 
data  below  250  m  altitude,  and  this  is  why  the  average  pro¬ 
file  is  somewhat  erratic  at  the  lowest  altitudes.  At  most 
of  the  other  altitudes  several  hundred  values  were  used  in 
the  average. 

Standard  deviations  of  the  H2O  data  were  also  calcu¬ 
lated  for  each  altitude.  Water  vapor  profiles  representing 
the  average  value  plus  one  standard  deviation  and  the  aver¬ 
age  minus  one  standard  deviation  for  each  altitude  are  also 
shown  in  Figure  19. 

The  water  vapor  concentrations  are  represented  in 
Figure  19  in  parts  per  10,000.  Table  6  gives  the  equivalent 
of  100  parts  per  10,000  in  other  units  which  are  commonly 
used  for  water  vapor. 

For  purposes  of  comparison,  the  water  vapor  profile  for 
the  midlatitude  summer  standard  atmosphere  [2]  is  plotted  in 
the  same  units  in  Figure  20.  Above  about  600  m  the  standard 
atmosphere  and  the  average  profile  from  the  data  are  roughly 
equivalent.  Below  600  m,  the  average  profile  is  somewhat 
drier  than  the  standard  atmosphere.  Also  shown  in  Figure  20 
is  the  midlatitude  summer  profile  shifted  to  result  in  100% 
relative  humidity  at  the  surface.  The  entire  profile  was 
shifted  by  the  amount  of  water  vapor  (parts  per  10,000) 
required  to  give  saturation  at  the  surface.  This  did  not 


36 


Average  of  oil  doto 


H20  parts  per  10»»4 


FIGURE  19.  WATER  VAPOR  PROFILE  OBTAINED  FP.OM  AVER¬ 
AGING  ALL  DATA  IN  FLIGHTS  1  THROUGH 
10,  AND  PROFILES  RESULTING  FROM  ADDING 
AND  SUBTRACTING  THE  STANDARD  DEVIATION 
OF  THE  DATA  AT  EACH  ALTITUDE. 
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TABLE  6.  CONVERSION  FACTORS  FOR  WATER  VAPOR  UNITS 


PARTS  PER  10,000 

g/m^ 

TORR 

Sea  Level 

100 

7 . 56 

7.60 

1  km  Altitude 

100 

6.73 

6.77 

1 

1 
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Widlotitude  Summer 


H20  parts  per  10**4 


FIGURE  20.  MIDLATITUDE  SUMMER  STANDARD  ATMOS- 
PHERE  WATER  VAPOR  PROFILE  (MLS) 

AND  THE  MIDLATITUDE  SUMMER  PROFILE 
SHIFTED  TO  GIVE  100%  RELATIVE 
HUMIDITY  AT  THE  SURFACE. 
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result  in  saturation  at  any  of  the  higher  altitudes  . 
although  relative  humidity  is  high  at  all  altitudes.  The 
standard  atmosphere  profile  for  100%  humidity  at  the  surface 
contains  considerably  more  water  vapor  than  the  average- 
plus-one-standard-deviation  profile  obtained  from  the  data. 


3.2  TOTAL  COLUMN  WATER  AMOUNTS 

The  total  amount  of  water  in  each  of  the  vertical 
columns  sampled,  the  precipitable  water,  can  be  calculated 
from  the  mixing  ratio  profiles.  The  precipitable  water,  h 
in  cm,  in  a  column  of  water  is  given  as 


h 


qy(  z) 


(1) 


where  is  the  surface  pressure,  H  is  the  pressure  scale 

height  g  is  the  acceleration  of  gravity  q^(z)  is  the  vol¬ 
ume  mixing  ratio  at  the  height  z  and  PH2O  is  the  density  of 
liquid  water.  Figures  21  to  31  display  the  precipitable 
water  values  calculated  using  (1)  for  each  of  the  fliahts- 
Flight  1  had  a  separate  downward  and  upward  leg  and  the 
results  from  the  legs  have  been  plotted  separately-  Also 
results  from  Flight  11  have  not  been  included.  The  reason 
is  that  Flight  11,  apparently,  corresponds  to  the  flioht 
back  to  NASA-Langley  and  is  not  in  the  area  sampled  during 
the  other  flights.  In  evaluating  the  results  in  Figures  21 
to  31  one  must  remember  that  the  vertical  columns  sampled 
varied  from  one  location  to  another.  Therefore,  the  varia¬ 
tion  seen  in  the  figures  is  not  entirely  a  measure  of 
changes  in  the  atmospheric  content  of  water  vapor  but  in¬ 
stead  represents  the  variation  in  the  depth  of  the  atmo¬ 
sphere  sampled. 
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Flight  1  Up 


Latitude 


FIGURE  21.  PRECIPITABLE  WATER  AMOUNTS  IN  THE 
VERTICAL  AIR  COLUMNS  SAMPLED  BY 
THE  NASA-DIAL  INSTRUMENT. 
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FIGURE  22.  PRECIPITABLE  WATER  AMOUNTS  IN  THE 
VERTICAL  AIR  COLUMNS  SAMPLED  BY 
THE  NASA-DIAL  INSTRUMENT. 
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Precipitable  Water  (cm) 
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FIGURE  23.  PRECIPITABLE  WATER  AMOUNTS  IN  THE 
VERTICAL  AIR  COLUMNS  SAMPLED  BY 
THE  NASA-DIAL  INSTRUMENT. 
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FIGURE  24.  PRECIPITABLE  WATER  AMOUNTS  IN  THE 
VERTICAL  AIR  COLUMNS  SAMPLED  BY 
THE  NASA-DIAL  INSTRUMENT. 
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FIGURE  25.  PPECIPITABLE  WATER  AMOUNTS  IN  THE 
VERTICAL  AIR  COLUMNS  SAMPLED  BY 
THE  NASA-DIAL  INSTRUMENT. 
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FIGURE  26.  PRECIPITABLE  WATER  AMOUNTS  IN  THE 
VERTICAL  AIR  COLUMNS  SAMPLED  BY 
THE  NASA-DIAL  INSTRUMENT. 
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FIGURE  27.  PRECIPITABLE  WATER  AMOUNTS  IN  THE 
VERTICAL  AIR  COLUMNS  SAMPLED  BY 
THE  NASA-DIAL  INSTRUMENT. 
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FIGURE  28.  PRECIPITABLE  WATER  AMOUNTS  IN  THE 
VERTICAL  AIR  COLUMNS  SAMPLED  IN 
THE  NASA-DIAL  INSTRUMENT. 
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FIGURE  29.  7RECIPITABLE  WATER  AMOUNTS  IN  THE 
VERTICAL  AIR  COLUMNS  SAMPLED  IN 
THE  NASA-DIAL  INSTRUMENT. 
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FIGURE  30.  PRECIPITABLE  WATER  AMOUNTS  IN  THE 
VERTICAL  AIR  COLUMNS  SAMPLED  IN 
THE  NASA-DIAL  INSTRUMENT. 
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FIGURE  31.  PRECIPITABLE  WATER  AMOUNTS  IN  THE 
VERTICAL  AIR  COLUMNS  SAMPLED  IN 
THE  NASA-DIAL  INSTRUMENT. 
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One  can  now  ask  the  question  if  these  results  can  be 
to  infer  information  about  the  atmosphere  that  can  be 
in  atmospheric  propagation  calculations. 

FITTING  THE  WATER  VAPOR  DATA 

Even  though  the  water  vapor  profiles  in  Figures  3  to  13 
a  fair  amount  of  structure  they  can  be  fit  reasonably 
with  a  relationship  of  the  form 

qy(  z)  =  Qq  +  m  logj^QZ  (  2  ) 

Table  7  gives  the  calculated  values  of  m  and  q^  for  the 
flights.  Flights  2  and  6  were  excluded  from  this  analysis 
due  to  the  dramatic  structure  (see  Figures  4  and  8).  The 
average  value  of  -31.005  is  similar  to  the  value  of  -29.8 
from  Starr  et  al  [5]  for  annually  and  latitud inally  averaged 
conditions  for  35°  N  during  the  IGY.  With  the  average  value 
of  m  one  can  solve  for  the  surface  values  Of  the  volume 
mixing  ratio,  qyo '  required  to  reproduce  the  average  pro¬ 
files  given  in  Figures  3  to  13.  Then,  with  the  surface 
value  known  one  can  solve  for  the  surface  temperature  re¬ 
quired  to  give  the  calculated  q^o*  These  surface  tempera¬ 
tures  can  then  be  compared  against  the  data  given  in  Tables 
2  to  4 . 

The  surface  value  of  q^^  is  given,  using  the  Claussius- 
Clapeyron  equation,  as 


used 

used 

3.3 

show 

well 


q 


VO 


r^  — exp  [5412 
o 


1 

273 


(3) 


where  r^  is  the  surface  value  of  the  relative 
is  the  surface  pressure  and  T^  is  the  surface 
Solving  for  Tq  gives 


humidity, 
temperature . 
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TABLE  7.  CALCULATED  SLOPE,  M,  AND  INTERCEPTS,  00,  TO  THE 
EMPIRICAL  FITS  TO  THE  AVERAGED  WATER  VAPOR 

PROFILES. 


FLIGHT 

SLOPE 

M 

INTERCEPTS 
(PP  10,000) 

1 

-43.185 

299 

3 

-19.932 

198 

4 

-23.253 

192 

5 

-13.288 

178 

7 

-23.253 

209 

8 

-43 . 185 

250 

9 

-33.219 

223 

10 

-43.185 

256 

11 

-36.541 

237 

I 

All  Flights 

-31 .005 

1 

227  ! 

. J 
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T  . 
o 


In  [ 


(4) 


5412 


6.11  r 

c 

P  q 
o  ^vo 


exp 


5412 

273 


)] 


The  values  of  and  calculated  from  (3)  and  (4) 

for  the  flights  are  given  in  Table  8.  In  the  calculations 
the  surface  relative  humidity  is  assumed  to  be  63%,  based  on 
the  results  in  Table  2.  The  calculated  surface  temperatures 
agree  well  with  those  measured  (see  Tables  2  to  4).  Using 
(3)  and  (4)  one  could  use  the  DIAL  data  to  infer  the  surface 
temperature  at  the  base  of  the  vertical  columns  sampled. 

Using  the  empirical  relationship  (2)  and  the  average 
slope  and  intercept  over  all  the  flights  the  surface  tem¬ 
peratures  required  to  produce  the  precipitable  water  in  the 
sampled  columns  of  air  were  iteratively  solved  for  using  the 
secant  method  [7].  Then,  with  the  surface  temperature  known, 
the  precipitable  water  in  a  vertical  column  extending  from 
the  surface  to  1361  m  was  calculated  as  well  as  the  average 
surface  temperature  for  each  flight.  The  results  are  given 
in  Tables  9  to  18.  The  surface  temperatures  inverted  from 
this  method  are  a  bit  higher  than  those  in  Tables  2  to  4  but 
are  in  reasonable  agreement.  The  precipitable  water  values 
are  also  consistent  with  climatological  values  for  that 
latitude  and  temperature  [5].  This  tends  to  indicate  that 
the  data  obtained  by  the  NASA  DIAL  instrument  were  similar 
to  climatological  conditions.  To  test  this  the  inversion  of 
surface  temperatures  was  redone  using  a  climatological  ex¬ 
pression  for  tlie  relative  humdity  suggested  by  Manabe  and 
Wetherald  [8] , 
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TABLE  8.  CALCULATED  VALUES  OF  THE  SURFACE 
VOLUME  MIXING  RATIO  AND  SURFACE 
FROM  THE  FLIGHT  AVERAGED  DATA. 


FLIGHT 

(pp  10,000) 

T 

(K) 

1 

299 

304 

3 

198 

298 

4 

192 

297 

5 

178 

296 

7 

209 

298 

8 

250 

301 

9 

223 

300 

10 

256 

302 

11 

.  237 

301 

All  Flights 

227 

300 

WATER  VAPOR 
TEMPERATURE 
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TABLE  9.  INVERTED  TEMPERATURES  AND  CALCULATED  TOTAL  COLUMN  PRECIPITABLE 
WATER  AMOUNTS  FOR  FLIGHT  1. 
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THE  AVERAGE  SURFACE  TEMPERATURE  IS  300.26 


TABLE  9.  (continued) 
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INVERTED  TEMPERATURES  AND  CALCULATED  TOTAL  COLUMN  PRECIPITABLE 
WATER  AMOUNTS  FOR  FLIGHT  2 . 


CNCNOcNinm^neo^OOncDD'^r^oeocD^rinf^ai'^oinoino^'^OoDOOCNOCNeo^r 

’•-CN’^'^CNTiinn^inTrin^^p^co^CNOOm^cocNOCN'^cN'^CNOtPiniDOCNCNOcxOtP 


CNCNCNCNCNCN*-OtNCN(N(NCNCN(N(NCNCNCNCN^CN*^?N(NCNCNCNCN(N^CNCNCN^CNCNCNCNCN^ 


OnoDO^CNcronTi<»-ncNinm^co^'^'^^nn^o^n^r)Tjr*»eonif)n'«T^CNTj'^tp 

■^ino)’^;pif)nr)cconr-TjTjcNa)if)cooor5CNOODmoDCNODCNeocN'^r-^OeoinoinooD 


oocoooo^oDQOQQO’^’^'^Oococo^cor^ocooooooeo'^Q'^ODaiooDOeocc 

ooooooooDOOOOOOCOOoooooooooooooooooooooooo 

CNCNCNCNCNnCNCNnnnDnnnnrOCNCNCNCNCNCNCNCNCNCNCNCNCNCNnnnCNCNCNCNCNCNCN 


^'*-tPTio)r^^icoD'*-mCT)cpcc^'^0^’^’^'^fNoo'^0’'TCN’«TrMTTor-ocomTjr'con-^o 

coooDCDODOccOCTOCTJOOOO'^^ronronncNnnnnnnncNOOOr^nninno’^ 

000000000^000 - '-000000000000000^0-0000000 


—  —  —  —  CCVTCOCCCCtT' 

U)U)lPiriPU)irU)<DU)U)U)lClCtDtCn^^nnTTTTTTTTTTTJTTTTTTTTCDtPlO<J1^TJTrTJTJTJ 

nronnr)r)nr)nr)r)r)nr5ror)r*'r-r'r^r'h'r'r^h'r*»r-r-r-r>“r*»r3r)r)r)r-r-r*-r*-r-r- 


IP 

IP 

IP 

IP 

IP 

IP 

IP 

IP 

oc 

00 

00 

n 

c 

n 

n 

tr 

n 

00 

iniomininmmmmmminmmmmminmininntTtT^tT^tT^^tTTynnnnnnnnr) 

CNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCN 

r^r^r^r^r^r^r*»r^r^r^r^r^r*»r^r^r^t'-.t'-.t'-.r-r-r'-.r*'r>‘r^r-r^r^r-r.r^r^r^r-r^r-r^r^r^r^r^ 


icinTjcN'-oa)r-icmnr)'-oo)cocDinnr)CNOOCDr-in^r5CN'-oeor^ 
cNCNCNCNCNOi'-'^'-^'-  —  mininmininininn^n 


n  (N  —  o  00  r* 
n  n  n  tr  n  n 


(r(rcrcr(^(r(c(r(r(r(r(r(r(ri0i0i£)i0(0u>i0(£)ininminininininminminininminmmin 

DroDDDnDnDnDnDnnnnnnronnnDnnDnonnnDroDnnnDnn 


iripiDiccoiru>uPoirir(ru?toioiotoioy)tDccioy?toy?ictoicictDVDCDtPvoiPU5^rvDtPtotp 
—  nm  —  Din  —  nin  —  Din  —  nm—  nm  —  nm  —  nm  —  Din  —  nm  —  nin  —  nin  —  nin'-n 


CNCscNnnDTjTjnintnininininr'h*r^cococoDa)0)000  —  —  —  CNCNCNDDD^rn^mm 
CNCNCNCNfNfMCNfNCNfNfNfNCNfSfNfNfNfNCNDCNfNDCNOiCNCNfNfNfNCNCNCNCNfNCNfNCNCNCNfS 


50 
O 


THE  AVERAGE  SURFACE  TEMPERATURE  IS  299  28 


TABLE  11.  INVERTED  TEMPERATURES  AND  CALCULATED  TOTAL  COLUMN  PRECIPITABLE 
WATER  AMOUNTS  FOR  FLIGHT  3. 
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THE  AVERAGE  SURFACE  TEMPERATURE  IS  298.96 


TABLE  12.  INVERTED  TEMPERATURES  AND  CALCULATED  TOTAL  COLUMN  PRECIPITABLE 

v/ater  amounts  for  flight  4  . 


60 


TMf  AVFffAGE  MJOfACF  TEMTEffAfUOE  t  297  48 


TABLE  13.  INVERTED  TEMPERATURES  AND  CALCULATED  TOTAL  COLUMN  PRECIPITABLE 
WATER  AMOUNTS  FOR  FLIGHT  5. 
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THE  AVERAGE  SURFACE  TEMPERATURE  IS  299  31 


TABLE  14.  INVERTED  TEMPERATURES  AND  CALCULATED  TOTAL  COLUMN  PRECIPITABLE 
WATER  AMOUNTS  FOR  FLIGHT  6 . 
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THE  AVERAGE  SURFACE  TEMPERATURE  IS  299  76 


TABLE  15.  INVERTED  TEMPERATURES  AND  CALCULATED  TOTAL  COLUMN  PRECIPITABLE 
WATER  AMOUNTS  FOR  FLIGHT  7 . 
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THE  AVERAGE  SURFACE  TEMPERATURE  IS  299.42 


TABLE  16.  INVERTED  TEMPERATURES  AND  CALCULATED  TOTAL  COLUMN  PRECIPITABLE 
WATER  AMOUNTS  FOR  FLIGHT  8. 
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THE  AVERAGE  SURFACE  TEMPERATURE  15  297  95 


TABLE  17.  INVERTED  TEMPERATURES  AND  CALCULATED  TOTAL  COLUMN  PRECIPITABLE 

water  amounts  for  flight  9 . 
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r(z)  =  -  0.02  )  /  0.98 


(5) 


o 


where  is  the  relative  humidity  at  the  surface  and  Pq  is 

the  surface  pressure.  With  (5)  the  mixing  ratio  is  given  by 


where  the  temperature  T(z)  is  given  by 


(7) 


=  Tq  -  r^z 


where  is  the  dry  adiabatic  lapse  rate.  Using  the  above 

approach  the  resultant  surface  temperatures  are  quite  simi¬ 
lar  to  those  obtained  using  the  empirical  relationship  de¬ 
rived  from  the  NASA  DIAL  data,  the  absolute  difference  being 
on  the  order  of  a  degree  K  or  less.  The  agreement  is  not 
unexpected  seeing  that  the  Manabe-Wetherald  derived  water 
vapor  profile  is  similar  to  that  obtained  from  the  midlati¬ 
tude  summer  atmosphere.  What  this  indicates,  then,  is  that 
one  could  use  the  DIAL  derived  precipitable  water  amounts 
and  an  assumed  profile  shape  to  invert  surface  temperatures. 
Conversely,  if  one  only  has  the  surface  temperature  and 
surface  relative  humidity  one  could  calculate  the  precipi¬ 
table  water  amount  reasonably  well  by  using  a  climato- 
logically  derived  water  vapor  profile  shape. 

The  inverted  absolute  temperatures  are  clearly  depen¬ 
dent  upon  the  assumptions  made  concerning  the  humidity  dis¬ 
tribution.  The  temperature  differences,  however,  are  less 
so  and  are  probably  more  useful  in  that  they  can  be  used  to 
locate  sea  surface  temperature  variations,  such  as  might  be 
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associated  with  the  Gulf  Stream.  The  temperatures  given  in 
Tables  9  to  18  show  temperature  variations  on  both  spatial 
and  temporal  scales  that  can  be  interpreted  as  relating  to 
water  temperature  variations. 

The  flight  pattern  used  during  the  tests  was  such  that 
the  plane  flew  the  same  flight  path  at  different  times  of 
the  day.  For  example,  the  downward  leg  of  flight  1  and 
flights  4  and  8  covered,  essentially,  the  same  flight  paths. 
The  inverted  temperatures  from  the  overlapping  portions  show 
a  variation  from  one  time  to  another  on  the  order  of  several 
degrees.  The  temperatures  inverted  from  flight  1  are  higher 
than  those  from  the  later  flights.  Likewise  the  tempera¬ 
tures  inverted  from  the  upward  leg  of  flight  1  are  higher 
than  those  frcxn  the  overlapping  portion  of  flight  5.  An 
oceanographic  analysis  provided  by  the  National  Weather 
Service  and  National  Earth  Satellite  Service  indicated  that 
the  Gulf  Stream  was  observed  to  be  in  the  vicinity  of  the 
test  area  and  the  variations  in  the  inverted  surface  temper¬ 
ature  could  be  the  result  of  Gulf  Stream  meanderings. 

The  surface  temperature  inversion  method  presented  here 
is  much  simpler  than  that  of  Rosenberg  and  Hogan  [9] .  Their 
approach  involved  using  a  three  frequency  DIAL  to  invert  the 
temperature  and  humidity  profiles.  In  their  feasibility 
study  they  used  H2O  lines  and  considered  both  horizontal  and 
vertical  viewing  angles.  Their  study  indicated  that  one  can 
invert  the  temperature  and  humidity  profiles  although  the 
error  in  temperature  and  humidity  increases  with  increasing 
horizontal  and  vertical  range.  Tiie  simple  approach  pre¬ 
sented  here  assumed  knowledge  about  the  shape  of  both  the 
temperature  and  humidity  profiles.  An  extensive  analysis  of 
the  sensitivity  of  the  inversion  method  to  the  model  assump¬ 
tions  has  not  been  performed  at  this  time. 
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4 


INFLUENCE  OF  H2O  PROFILES  DERIVED  FROM  THE  DIAL  DATA 
ON  DF  LASER  PROPAGATION 


4.1  MODEL  FOR  CALCULATION  OF  DF  LASER  TRANSMITTANCE 

4.1.1  TRANSMITTANCE  FOR  A  LASER  POWER  SPECTRUM 

An  average  transmittance  over  a  laser  power  spectrum 
can  be  defined  in  terms  of  the  transmittances  at  the  indi¬ 
vidual  laser  frequencies.  The  transmittance  at  some  range  r 
is  the  ratio  of  the  total  beam  power  P  at  that  range  to  the 
beam  power  Pg  at  the  source: 


T  =  P/Pg 


The  total  beam  power  at  range  r  is  the  sum  of  the  source 
powers  in  the  individual  lines  times  the  transmittances  for 
the  individual  lines: 


T 


If.  P 
11s 


T. 

1 


where  the  f^^  define  the  power  spectrum  of  the  laser  and  T^^ 
is  the  transmittance  for  the  ith  line.  The  laser  transmit¬ 
tance  is  therefore  the  sum  of  the  individual  transmittances 
weighted  by  the  power  spectrum: 
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4.1.2  TRANSMITTANCE  CALCULATION  FOR  THE  LIDAR  DATA 

The  lidar  data  consist  of  water  vapor  amounts  in  parts 
per  ten-thousand  at  each  of  various  altitudes.  Typically, 
there  were  data  for  60  layers  between  200  m  and  1400  m  alti¬ 
tude.  Transmi ttances  were  calculated  for  these  water  vapor 
profiles,  assuming  midlatitude  summer  standard  atmosphere 
conditions  for  everything  except  water  vapor. 


The  transmittance  calculations  were  done  both  for  a  cw 
and  for  a  pulsed  DF  laser.  For  the  cw  laser  calculations 

the  power  spectrum  of  the  MIRACL  laser  was  used,  and  it  is 

assumed  that  this  is  representative  of  high-powered  cw  DF 
lasers.  This  power  spectrum  is  given  in  Table  19  [10].  For 
the  pulsed  laser  calculations,  the  power  spectrum  listed  in 
Table  20  was  used.  This  is  not  a  power  spectrum  for  an 

actual  device,  but  has  the  general  characteristics  (more 
power  at  higher  v  and  J  levels)  of  high-power  pulsed  DF 

lasers.  In  the  pulsed  power  spectrum,  it  is  assumed  that 
the  lines  which  occur  in  the  region  of  the  4.3  m  CO2  band 
are  removed,  but  that  lines  at  longer  wavelengths  are  still 
present.  This  is  important,  because  these  longer-wavelength 
lines  suffer  some  absorption  by  H2O.  This  is  not  an  issue 
for  cw  lasers,  since  they  operate  on  the  shorter-wavelength 
transitions  and  transmit  almost  no  energy  in  the  region  of 
the  CO2  band  and  beyond . 

A  computer  code  named  TRANS  was  written  to  perform  DF 
laser  transmittance  calculations.  TRANS  was  used  to  calcu¬ 
late  transmittance  at  each  laser  frequency  for  each  layer, 
and  these  results  were  combined  to  produce  transmi ttances 
for  2-km  paths  oriented  vertically  and  at  BO-degree  and  60- 
degree  zenith  angles.  The  data  for  the  lowest  layer  were 
also  used  to  estimate  horizontal  transmittance  at  the  sur¬ 
face  for  a  2-km  range. 
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table  19  .  CW  DF  LASER  POWER  SPECTRUM  USED  IN 
TRANSMITTANCE  CALCULATIONS. 


LINE  POWER  SPECTRUM 


Pi  ( 8  ) 

0.07 

Pl(9) 

0 . 31 

Pi ( 1 0 ) 

0.16 

Pi (11  ) 

0.04 

P2  (  8  ) 

0.11 

^2(9) 

0.08 

P2(10) 

0.07 

P2(ll) 

0.01 

P2(12) 

0.01 

P3(7) 

0.03 

P3  (  8  ) 

0.04 

P3(9) 

0 . 06 

P3(10) 

0.01 
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TABLE  20 


PULSED  DF  LASER  POWER  SPECTRUM 
USED  IN  TRANSMITTANCE  CALCULATIONS 


LINE 

POWER  SPECTRUM 

Pi  (8) 

0.033 

Pl(9) 

0.041 

Pl(lO) 

0.033 

P2(6) 

0.033 

P2(8) 

0.066 

P2(9) 

0.049 

P2(10) 

0.033 

P3(6) 

0.041 

P3(7) 

0.074 

P3(8) 

0.082 

P3(9) 

0.074 

P4{6) 

0.049 

P4(7) 

0.066 

P4(8) 

0.049 

P4  ( 9 ) 

0.033 

P5(6) 

0.041 

P5(7) 

0.057 

P8(6) 

0 . 049 

0.025 

P9(5) 

0.016 

P9(6) 

0.041 

P9(7) 

0.016 
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Each  calculation  by  TRANS  of  transmittance  at  a  laser 
frequency  in  a  layer  of  atmosphere  consisted  of  four  parts: 

-  line  absorption  by  H2O; 

-  line  absorption  by  atmospheric  species  other  than 

H2O; 

-  H2O  continuum;  and 

-  N2O  continuum. 

Each  of  these  is  discussed  below. 

4. 1.2.1  Line  Absorption  By  H^O 

Volume  absorption  coefficients  for  each  laser  freauency 
in  the  power  spectrum  for  midlatitude  summer  conditions  at 
sea  level  and  at  5  km  were  obtained  from  the  tabulation  by 
Manning  and  Matise  [11]  .  These  coefficients  were  converted 
to  cross  sections  ( cm'^ /molecule )  ,  which  are  independent  of 
water-vapor  amount.  The  actual  cross  section  at  the  altitude 
of  each  individual  layer  was  found  by  logarithmic  interpola¬ 
tion  of  the  values  at  0  and  5  km.  The  optical  depth  in  any 
layer  is  the  product  of  the  absorption  cross  section,  the 
number  density  of  H2O  molecules  in  the  layer,  and  the  path 
length  through  the  layer. 

4. 1.2. 2  Line  Absorption  By  Species  Other  Than  H2O 

Volume  absorption  coefficients  for  species  other  than 
H2O  were  obtained  from  the  tabulations  of  Manning  and  Matis° 
[11]  for  a  midlatitude  summer  standard  atmosphere  at  alti¬ 
tudes  of  0  and  5  km.  The  absorption  coefficients  tabulated 
by  Manning  and  Matise  were  based  on  line-by-line  calcula¬ 
tions  which  included  a  sub-Lorentz  lineshape  for  CO2  and 
far-wing  effects  for  all  major  species. 
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The  actual  absorption  coefficient  in  any  level  is  found 
by  logarithmic  interpolation  between  the  values  at  0  and  5 
km.  The  optical  depth  of  the  layer  for  non-H20  line  ab¬ 
sorption  is  the  product  of  the  volume  absorption  coefficient 
and  the  path  length  through  the  layer. 

4. 1.2. 3  Water-Vapor  Continuum 

The  optical  depth  in  a  layer  due  to  the  water-vapor 
continuum  is  calculated  using  an  algorithm  based  on  the 
measurements  of  White,  et  al.  [12]  and  Watkins,  et  al.  [13]. 
The  algorithm  calculates  absorption  as  a  function  of  fre¬ 
quency,  temperature,  and  H2O  concentration.  It  is  described 
in  the  High  Energy  Laser  Propagation  Handbook  [10]  . 

4. 1.2. 4  Nitrogen  Continuum 

The  optical  depth  through  the  N2  continuum  is  calcu¬ 
lated  for  each  layer  and  for  each  laser  frequency  using  an 
algorithm  based  on  the  data  of  Shapiro  and  Gush  [14].  The 
algorithm  calculates  absorption  by  the  N2  continuum  as  a 
function  of  frequency,  temperature,  and  atmospheric  density. 
It  is  described  in  the  High  Energy  Laser  Propagation  Hand¬ 
book  [10]. 

4.2  RESULTS  OF  TRANSMITTANCE  CALCULATIONS 

Computer  program  TRANS,  described  in  Section  4.1,  was 
used  to  calculate  transmi ttances  for  each  measured  water 
vapor  profile.  The  data  did  not  extend  to  the  surface.  In 
order  to  do  these  calculations,  the  water  vapor  concentra¬ 
tion  at  the  lowest  altitude  in  each  profile  was  extended 
down  to  the  surface,  and  the  concentration  in  the  highest 
altitude  was  extended  up  to  2  km.  A  more  realistic  method 
of  extending  the  measured  profiles  in  both  directions  could 
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certainly  have  been  found;  in  fact,  for  these  data,  the 
first  and  last  points  in  each  profile  are  probably  the  least 
accurate.  The  laser  transmittances  were  so  high,  however, 
that  it  was  decided  that  such  "fine  tuning"  of  the  model  was 
not  necessary. 

The  transmittance  calculations  were  done  for  two-km 
ranges  for  vertical  paths,  paths  at  30-degree  and  60-degree 
zenith  angles,  and  for  horizontal  paths.  Since  the  same 
range  was  used  for  each  path,  it  was  possible  for  transmit¬ 
tance  to  actually  be  higher  for  the  slant  paths  than  for  the 
vertical  path.  This  would  occur  if  there  was  a  high  concen¬ 
tration  of  water  vapor  in  the  upper  layer;  the  vertical  path 
would  pass  through  that  water  vapor,  but  the  slant  paths 
would  have  ended  before  reaching  the  upper  layers.  This 
actually  did  occur  several  times. 

The  sea-level  horizontal  path  transmittance  calcula¬ 
tions  were  generally  not  meaningful.  An  example  is  shown  in 
Figure  32  for  flight  5.  These  calculations  were  based  on 
the  data  point  for  the  lowest  altitude  in  each  profile;  as 
just  noted,  this  point  is  not  expected  to  be  very  accurate. 
The  strong,  random  fluctuations  in  Figure  32  indicate  a  high 
noise  level  in  the  data  for  the  lowest  altitude  of  each 
profile. 

4.2.1  CW  LASER  TRANSMITTANCE  FOR  INDIVIDUAL  PROFILES 

The  results  of  the  cw  laser  transmittance  calculations 
for  flight  1  are  shown  in  Figure  33  for  the  vertical  path 
and  the  two  slant  paths.  There  are  two  curves  on  each  plot 
because  the  aircraft  changed  directions  during  the  course  of 
flight  1.  The  directions  are  shown  by  arrows  on  the  first 
plot  of  the  figure.  The  transmittances  are  lower  for  the 
30-degree  and  60-degree  zenith  angle  paths.  This  is 
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FJ»gM  5  9/24/B2  horizontoi  poth 


lotttLriM 


FIGURE  32,  CW  DF  LASER  TRANSMITTANCE  FOR  A 
HORIZONTAL  PATH  AT  SEA  LEVEL 
CALCULATED  FROM  THE  DATA  FOR 
FLIGHT  5. 
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y  9/24/82  v*ft>coi  poif* 


FIGURE  33.  CW  DF  LASER  TRANSMITTANCE  CALCULATED 
FROM  DATA  FOR  FLIGHT  1. 
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generally  the  expected  result,  despite  the  fact  that  all 
three  paths  have  the  same  range;  the  water  vapor 
concentration  generally  falls  off  with  altitude,  and  the 
slant  paths  have  longer  paths  through  the  lower  layers  where 
the  concentration  is  high.  The  difference  between  the 
vertical  transmittance  and  the  transmittance  for  the  60- 
degree  zenith  angle  ranges  from  zero  to  about  4% 
(transmittance  units). 


The  calculation  results  for  cw  lasers  for  flight  2  are 
shown  in  Figure  34.  In  some  regions  the  transmittances  are 
nearly  the  same  for  the  vertical  and  slant  paths,  while  for 
other  regions  the  transmittance  for  the  60-degree  zenith 
angle  is  as  much  as  5%  lower  (transmittance  units).  Trans¬ 
mittance  for  2-km  paths  is  around  90%. 


Cw  results  for  flight  3  are  shown  in  Figure  35  for  the 
vertical  path.  The  slant-path  transmittances  (not  shown) 
were  lower,  but  only  by  0  to  2%.  The  same  was  true  for 
flights  4  and  5  (Figures  36  and  37).  The  slant  path  trans¬ 
mittances  for  flights  6  and  7  (Figures  38  and  39)  were 

nearly  identical  to  the  vertical  transmittances.  For  flight 
8  (Figure  40)  the  slant  path  transmittances  were  slightly 

lower  in  some  places  and  about  the  same  as  the  vertical 

transmittances  in  other  places.  Some  of  the  slant  path 
transmittances  for  flight  9  (Figure  41)  were  higher  than  the 
corresponding  vertical  path  transmittances,  indicating  high 
water  vapor  concentrations  in  the  upper  layers.  For  flights 
10  and  11  (Figures  42  and  43)  the  slant  path  transmittances 
were  slightly  lower  than  the  vertical  path  values.  The 

aircraft  changed  directions  during  flight  11,  causing  the 
loop  in  the  transmi ttance-versus-longitude  plots  of 
Figure  43  . 
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FIGURE  34.  CW  DF  LASER  TRANSMITTANCE  CALCULATED  FROM 
DATA  FOR  FLIGHT  2. 
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Fkgnt  J  (/24/S2  vtrticoi  path 


FIGURE  35.  CW  DF  LASER  TRANSMITTANCE  CALCULATED 
FROM  DATA  FOR  FLIGHT  3. 


Flight  4  9/24/B2  vtrticol  poth 


FIGURE  36.  CW  DF  LASER  TRANSMITTANCE  CALCULATED 
FROM  DATA  FOR  FLIGHT  4. 
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Fli9^t  b  9/24/62  vtrticol  poth 


FIGURE  37.  CW  DF  LASER  TRANSMITTANCE  CALCULATED 
FROM  DATA  FOR  FLIGHT  5. 


FJight  6  9/24/62  v«ftico*  poth 


FIGURE  38.  CW  DF  LASER  TRANSMITTANCE  CALCULATED 
FROM  DATA  FOR  FLIGHT  6. 
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Fk«nt  7  •/24/B2  wbcol  »ath 


FIGURE  39 .  CW  DF  LASER  TRANSMITTANCE  CALCULATED 
FROM  DATA  FOR  FLIGHT  7. 
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FIGURE  40.  CW  DF  LASER  TRANSMITTANCE  CALCULATED 
FROM  DATA  FOR  FLIGHT  8. 
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FIGURE  41.  CW  DF  LASER  TRANSMITTANCE  CALCULATED 
FROM  DATA  FOR  FLIGHT  9. 
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FIGURE  42.  CW  DF  LASER  TRANSMITTANCE  CALCULATED 
FROM  DATA  FOR  FLIGHT  10. 
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FIGURE  43.  CW  DF  LASER  TRANSMITTANCE 
CALCULATED  FROM  DATA  FOR 
FLIGHT  11. 
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The  transmittances  plotted  in  Figures  33  through  43  are 
all  around  90%  for  a  two-km  range.  These  calculations  in¬ 
cluded  line  and  continuum  absorption  by  all  atmospheric 
gases.  Although  there  is  some  fluctuation  in  the  transmit¬ 
tance,  which  must  be  caused  by  water  vapor  variation,  the 
fluctuations  are  small,  and  no  spatial  pattern  is  seen  in 
the  results.  The  gulf  stream  boundary  appears  to  have  no 
influence  on  cw  DF  laser  transmittance. 


4.2.2  PULSED  LASER  TRANSMITTANCES  FOR 
PROFILES 


INDIVIDUAL 


Transmittance  calculations  were  also  done  for  a  pulsed 
DF  laser  for  each  measured  profile.  The  transmittances  were 
somewhat  lower  than  the  cw  transmittances.  Part  of  the  rea¬ 
son  for  this  is  that  several  laser  lines  on  the  long-wave¬ 
length  side  of  the  4.3  ym  CO2  band  are  included  in  the 
pulsed  laser  power  spectrum,  and  these  lines  are  subject  to 
fairly  strong  H2O  absorption.  Whether  these  lines  would  be 
present  in  the  power  spectrum  of  an  actual  device  would  de¬ 
pend  on  how  the  device  was  configured,  and  in  particular  how 
the  lines  in  the  region  of  the  CO2  band  were  suppressed. 


The  results  of  the  pulsed  DF  laser  transmittance  calcu¬ 
lations  for  flight  1  are  shown  in  Figure  44.  As  was  noted 

earlier,  the  aircraft  changed  directions  during  flight  1,  so 
that  each  latitude  was  covered  twice.  The  slant-path  trans¬ 
mittances  are  slightly  lower  than  the  vertical-path  trans¬ 
mittances.  The  calculated  transmittances  for  a  pulsed  DF 
laser  for  flights  2  through  10  are  shown  in  Figures  45 

through  53.  In  each  of  these  figures,  the  transmittances 

for  the  vertical  path  and  for  the  30-degree  and  60-degree 
zenith  angles  are  shown  on  the  same  plot.  For  flight  2 

(Figure  45)  the  slant-path  transmittances  are  lower  than  the 
vertical-path  transmittances  for  part  of  the  flight  and  are 
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Fbght  1  9/24/62  v^rttcol  poth 


FKght  1  9/24/62  30  dmof  zsnith  ooo^ 


FkgAt  I  9/24/62  60  zsmtn  ongi* 


l«t I  tod* 

FIGURE  44.  PULSED  DF  LASER  TRANSMITTANCE 

CALCULATED  FROM  DATA  FOR  FLIGHT  1 
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Fllgri  2  9/24/62 


FIGURE  45.  PULSED  DF  LASER  TRANSMITTANCE 

CALCULATED  FROM  DATA  FOR  FLIGHT  2 


Flight  3  3/24/82 


FIGURE  46.  PULSED  DF  LASER  TRANSMITTANCE 

CALCULATED  FROM  DATA  FOR  FLIGHT  3 
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Fll^l  4  9/24/62 


FIGURE  47.  PULSED  DF  LASER  TRANSMITTANCE 

CALCULATED  FROM  DATA  FOR  FLIGHT  4 


FMgM  S  9/24/82 


FIGURE  48.  PULSED  DF  LASER  TRANSMITTANCE 

CALCULATED  FROM  DATA  FOR  FLIGHT  5 
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Fiiint  t  a/24/(2 


FIGURE  49.  PULSED  DF  LASER  TRANSMITTANCE 

CALCULATED  FROM  DATA  FOR  FLIGHT  6 


Flight  7  9/24/B2 


FIGURE  50.  PULSED  DF  LASER  TRANSMITTANCE 

CALCULATED  FROM  DATA  FOR  FLIGHT  7 
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Flight  ■  S/24/(2 


FIGURE  51.  PULSED  DF  LASER  TRANSMITTANCE 

CALCULATED  FROM  DATA  FOR  FLIGHT  8 


Flight  3  3/24/82 


FIGURE  52.  PULSED  DF  LASER  TRANSMITTANCE 

CALCULATED  FROM  DATA  FOR  FLIGHT  9 
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FIGURE  53.  PULSED  DF  LASER  TRANSMITTANCE 

CALCULATED  FROM  DATA  FOR  FLIGHT 

10  . 
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about  the  same  for  other  parts;  at  one  point  the  slant-path 
transmi ttances  are  higher.  As  noted  in  the  discussion  of 
the  cw  laser  transmittance  calculations,  this  depends  on  the 
shape  of  the  water  vapor  profile.  For  most  of  the  remaining 
flights,  the  slant-path  transmi ttances  were  lower  than  the 
vertical-path  transmittances ,  although  the  curves  do  cross 
briefly  for  flights  5  and  6. 


The  transmittances  calculated  for  a  pulsed  DF  laser 
were  generally  in  the  range  of  85%  to  88%  for  a  2-km  range. 
As  with  the  cw  transmittances,  the  pulsed  laser  transmit¬ 
tances  contain  small  fluctuations  which  must  be  due  to  vari¬ 
ations  in  the  water  vapor  profile  (because  other  parameters, 
such  as  CO2  concentration,  were  taken  from  standard  atmo¬ 
sphere  models  and  did  not  vary  from  sounding  to  sounding). 
There  was  no  significant  spatial  variation  in  the  results 
which  might  correlate  with  the  gulf  stream  boundary. 


4.2.3  TRANSMITTANCES  FOR  AVERAGE  AND  EXTREME  PROFILES 

Transmittances  were  calculated  using  TRANS  for  the 
average  profile  discussed  in  Section  3.  This  profile  was 
determined  by  taking  the  average  of  all  of  the  data  at  each 
altitude.  High  and  low  humidity  profiles  were  defined  by 
adding  and  subtracting  the  standard  deviation  of  the  data  at 
each  altitude  to  the  average  at  that  altitude.  Calculations 
were  also  done  for  these  profiles.  The  calculations  are  the 
same  as  those  done  for  the  individual  profiles:  transmit¬ 
tances  for  vertical  paths  and  30-degree  and  60-degree  zenith 
angles,  all  with  identical  ranges  of  2  km.  The  results  of 
these  calculations  are  shown  in  Table  21.  The  cw  transmit¬ 
tances  are  about  5%  higher  (in  transmittance  units)  than  the 
pulsed  transmittances.  The  difference  between  the  high- 
humidity  and  low-humidity  profile  transmittances  is  only  4% 
or  5%  transmittance. 
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TABLE  21 


DF  LASER  TRANSMITTANCES  CALCULATED  FOR  THE  AVERAGE 
WATER  VAPOR  PROFILE  DERIVED  FROM  THE  DATA,  AND  FOR 
PROFILES  DEFINED  BY  THE  AVERAGE  MINUS  AND  THE  AVERAGE 
PLUS  THE  STANDARD  DEVIATION  a  OF  THE  DATA  AT  EACH 
ALTITUDE. 


PROFILE 

AVE  -  0 

AVERAGE 

AVE  +  0 

CW  DF  LASER 

Vertical 

0 . 928 

0 . 910 

0 . 888 

30°  zenith 

0.925 

0.906 

0.884 

60°  zenith 

0.910 

0.890 

0 . 866 

PULSED  DF  LASER 

Vertical 

0 . 877 

0 . 860 

0 . 839 

30°  zenith 

0.873 

0.855 

0.834 

60°  zenith 

0.855 

0 . 836 

0 . 814 
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4.2.4  TRANSMITTANCE  CALCULATIONS  FOR  A  STANDARD  ATMO¬ 
SPHERE 

TRANS  was  used  to  calculate  cw  and  pulsed  DF  laser 
transmittances  for  2-km  paths  through  the  midlatitude  summer 
standard  atmosphere  and  through  the  midlatitude  summer  at¬ 
mosphere  with  water  vapor  shifted  so  that  the  surface  humid¬ 
ity  was  100%.  The  water  vapor  profile  for  this  modified 
midlatitude  summer  atmosphere  is  almost  identical  to  that  of 
the  tropical  standard  atmosphere. 

The  results  of  these  calculations  are  given  in 
Table  22.  The  transmittances  are  slightly  lower  than  those 
calculated  for  the  data,  because  the  water  vapor  amounts  are 
somewhat  higher  for  the  standard  atmosphere  than  for  the 
measured  profiles.  As  was  noted  for  the  calculations  for 
the  data  averages,  a  fairly  substantial  change  in  the  water 
vapor  amount  produces  only  a  small  change  in  the  transmit¬ 
tance  . 

4.3  DESCRIPTION  OF  THERMAL  BLOOMING  CALCULATIONS 

A  high  energy  laser  beam  which  propagates  through  the 
atmosphere  has  an  optimum  power  level  at  which  the  intensity 
at  the  target  is  maximized.  If  the  laser  power  is  increased 
above  the  optimum  value,  the  resulting  increase  in  thermal 
blooming  will  cause  a  decrease  in  the  intensity  at  the  tar¬ 
get.  The  optimum  power  depends  on  the  laser  power  spectrum, 
the  amount  of  atmospheric  absorption  on  each  laser  line,  and 
the  degree  of  focusing.  Optimum  power  for  a  given  laser 
system  and  scenario  can  be  calculated  using  the  BRLPRO  code. 

BRLPRO  is  a  scaling-law  nonlinear  propagation  code 
originally  developed  by  H.J.  Breaux  of  the  Army  Ballistic 
Research  Laboratory  [15]  .  It  is  in  common  use  in  the  prop- 
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TABLE  22 


.  DF  LASER  TRANSMITTANCES  CALCULATED  FOR  A  MIDLATITUDE 
SUMMER  STANDARD  ATMOSPHERE  AND  FOR  A  MIDLATITUDE  SUM¬ 
MER  ATMOSPHERE  WITH  THE  WATER  VAPOR  PROFILE  SHIFTED 
TO  RESULT  IN  100%  RELATIVE  HUMIDITY  AT  THE  SURFACE. 


PROFILE 

MIDLATITUDE 

SUMMER 

MIDLATITUDE  SUMMER 
WITH  100%  SURFACE 
HUMIDITY 

CW  DF  LASER 

Vertical 

0.898 

0 . 854 

30°  zenith 

0.892 

0.848 

60°  zenith 

0 . 876 

0.829 

PULSED  DF  LASER 

Vertical 

0.849 

0.806 

30°  zenith 

0.842 

0.799 

60°  zenith 

0 . 822 

0.778 
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agation  community,  and  is  now  incorporated  in  the  Armiy's 
EOSAEL  library  of  propagation  codes  [16].  The  code  performs 
calculations  for  a  beam  which  is  assumed  to  have  a  Gaussian 
cross  section.  It  steps  along  the  beam  path,  calculating 
the  beam  size  and  intensity  at  each  point  along  the  path, 
considering  atmospheric  absorption,  scattering,  turbulence, 
jitter,  and  thermal  blooming.  The  calculations  are  done  for 
a  beam  consisting  of  a  single  laser  line,  so  a  power  spec¬ 
trum  must  be  simulated  using  power-weighted  absorption  coef¬ 
ficients. 


The  calculations  described  here  were  done  for  a  high 
energy  laser  system  used  in  a  ship  defense  scenario.  It  is 
important  for  such  a  system  to  operate  near  its  optimum 
power,  and  ideally  the  concentrations  of  absorbing  atmos¬ 
pheric  species  should  be  known  at  all  altitudes  so  that  the 
optimum  power  level  can  be  determined.  In  the  absence  of 
such  atmospheric  profile  measurements,  standard  atmosphere 
profiles  would  probably  be  used  in  conjunction  with  surface 
measurements  to  determine  optimum  power.  The  calculations 
were  done  to  determine  the  difference  between  the  "real 
optimum  power"  based  on  the  measured  water  vapor  profiles 
determined  by  the  DIAL  system  and  the  "assumed  optimum 
power"  based  on  standard  atmosphere  information. 


A  schematic  of  the  scenario  is  shown  in  Figure  54.  A 
ship  is  the  target  of  an  incoming  missile.  A  DF  laser  sys¬ 
tem  is  on  an  escort  ship  located  1  km  from  the  targeted 
ship.  The  missile  is  currently  located  such  that  its  range 
from  the  laser  is  2  km  and  the  beam  path  is  elevated  from 
horizontal  by  30  degrees.  Both  cw  and  pulsed  DF  lasers  were 
considered.  The  source  aperture  diameter  is  1  m  and  the 
laser  is  capable  of  operating  at  powers  in  excess  of  1  Mw . 
For  the  pulsed  laser  calculations,  a  pulse  duration  of  2  ms 
and  a  repetition  frequency  of  50  hz  were  used.  Turbulence 
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was  included  in  the  calculations  with  of  1x10“^^  . 

A  5  m/s  cross  wind  was  assumed. 

Absorption  coefficients  for  both  the  measured  profiles 
and  the  standard  atmosphere  profiles  were  determined  from 
the  transmi ttances  for  cw  and  pulsed  systems  described  in 
Section  4.3  of  this  report.  These  absorption  coefficients 
are  power-weighted  values  averaged  over  the  beam  path. 

4.4  RESULTS  OF  THERMAL  BLOOMING  CALCULATIONS 

The  results  of  the  BRLPRO  calculations  are  shown  in 
Table  23  for  both  cw  and  pulsed  systems  operating  over  a  2- 
km  range.  The  calculations  were  done  for  the  average  pro¬ 
file  which  was  derived  from  the  DIAL  data  by  averaging  all 
of  the  data  at  each  altitude,  and  the  "ave+a"  and 
"ave- a"  profiles  which  represent  the  average  value  plus  and 
minus  the  standard  deviation  of  the  data  at  each  altitude. 
The  results  for  the  average  profile  show  the  typical  error 
which  might  have  been  expected  by  using  standard  atmospheres 
rather  than  actual  profiles  to  determine  optimum  laser  power 
on  the  day  the  data  were  taken.  The  results  for  the  "ave-a" 
and  "ave+a"  profiles  show  the  range  of  error.  The  "ave+a" 
profile  was  the  closest  of  the  three  to  the  standard  at¬ 
mosphere  . 

The  first  column  in  Table  23  contains  the  absorption 
coefficients  used  for  the  calculations.  As  noted  previous¬ 
ly,  these  were  obtained  from  the  transmittances  calculated 
earlier,  and  represent  power-weighted  coefficients  averaged 
over  the  path  for  each  power  spectrum  and  each  profile.  The 
second  column  lists  the  "assumed  optimum"  power  which  was 
obtained  through  a  BRLPRO  calculation  using  the  midlatitude 
summer  standard  atmosphere  profiles  for  all  species,  includ¬ 
ing  water  vapor.  The  third  column  contains  the  "real 
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TABLE  23.  INTENSITY  LOSS  CAUSED  BY  USING  STANDARD  ATMOSPHERE  PROFILES  RATHER  THAN  REAL  PROFILES 
DETERMINE  OPTIMUM  LASER  POWER. 
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optimum"  power  values  obtained  by  using  the  actual  water 
vapor  profile.  Note  that  the  pulsed  laser  can  operate  at 
much  higher  powers,  despite  the  slightly  higher  absorption 
coef f icients . 

Column  four  lists  the  target  intensities  which  would  be 
achieved  if  the  laser  operated  at  the  "assumed  optimum" 
power  in  an  atmosphere  with  the  appropriate  "real"  water 
vapor  profile.  This  is  the  intensity  which  might  be  ex¬ 
pected  if  no  equipment  were  available  to  measure  the  water 
vapor  profile  and  a  standard  atmosphere  were  assumed.  Column 
five  gives  the  target  intensities  obtained  using  the  "real 
optimum"  power;  this  is  the  maximum  target  intensity  which 
could  be  achieved  for  those  atmospheric  conditions.  The 
last  column  lists  the  percentage  loss  of  intensity  caused  by 
using  the  standard  atmosphere  rather  than  the  real  profile 
to  determine  optimum  power  for  the  laser. 

In  the  worst  case,  there  is  only  a  4%  loss  of  target 
intensity  caused  by  using  the  standard  atmosphere  rather 
than  the  real  water  vapor  profile  to  determine  the  optimum 
operating  power  for  the  laser.  Similar  calculations  were 
done  for  a  6-km  range,  with  similar  results. 
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5 

REDUCTION  AND  ANALYSIS  OF  THE  LFLID  SUPPORT  AND 
RADIOSONDE  TRANSMITTANCE  DATA 


Records  of  the  radiosonde  data  collected  during  the 
times  that  the  most  reliable  long-path  transmittance  mea¬ 
surements  were  made  (the  times  are  listed  in  Table  5, 
Section  2.2)  were  provided  to  OMI  by  John  Toner  of  the  Naval 
Air  Development  Center  (NADC),  the  organization  sponsoring 
the  LFLID  test  series.  Temperature  and  relative  humidity 
values  from  the  radiosonde  data  were  extracted  for  300  m 
intervals  between  sea  level  and  the  AMOS  observatory  alti¬ 
tude  of  3.05  km.  These  data  were  used  as  the  basis  of 
slant-path  transmittance  calculations  performed  using  the 
LOWTRAN  6  computer  model  [2]. 

Figures  55-72  present  plots  of  the  results  of  the 
LOWTRAN  calculations  together  with  the  correspondino  exper¬ 
imental  measurements.  The  figures  are  divided  into  two 
groups  corresponding  to  the  two  wavelength  regions  of  inter¬ 
est.  Figures  55-62  show  the  results  for  the  spectral  region 
between  3.0  and  6.0  urn.  Figures  63-67  show  the  comparisons 
between  the  calculated  and  measured  transmi ttances  for  the 
7.0  to  14.0  wm  region.* 

In  the  comparisons  shown  in  Figures  55-67  the  identifi¬ 
cation  code  of  the  measured  transmittance  spectrum  is  indi¬ 
cated  (see  Table  5  in  Section  2.2  for  an  explanation  of  this 
code).  The  LOWTRAN  calculations  shown  in  Figures  55-67  were 
performed  omitting  the  contribution  due  to  aerosol  atten- 


*  A  data  reduction  error  was  discovered  after  the  plots  were 
prepared.  Corrected  scale  values  for  the  7,0  to  14.0  r 
experimental  spectra  are  given  in  the  figure  captions. 
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trorwmfttonc* 


FIGURE 


FIGURE 


S06R1  and  LOWTRAN  (no  o«ro«ol) 


55.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 

SPECTRUM  S06R1  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION. 


S06R5  ond  LOWTRAN  (no  o«ro«ol) 


•ovalongth  (um) 


56.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 

SPECTRUM  S06R5  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION. 
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S12R1  and  LOWTRAN  (no  oorotol) 


FIGURE 


FIGURE 


57.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 

SPECTRUM  S12R1  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION. 


3.0  3.5  4.0  4.5  5.0  5.5  6.0 

wovelongtti  (unt) 


58.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 

SPECTRUM  S23R1  COMPARED  WITH  LOWTRAN  5  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION. 
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trontmUtonc* 


S23RJ  (MXi  LOWTRAN  (no  oorotol) 


FIGURE  59.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 

SPECTRUM  S23R3  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION. 


S2SR1  and  LOWTRAN  (no  oaroaol) 


wovatangth  (um) 


FIGURE  60.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 

SPECTRUM  S25R1  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION. 
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S25R4  and  LOWTRAN  (no  o*ro*ol) 


FIGURE 


FIGURE 


61  .  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 

SPECTRUM  S25R4  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION. 


S25R6  ood  LOWTRAN  (no  oorosol) 


o 


vrovelonyth  (um) 


62.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 

SPECTRUM  S25R6  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION. 
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Irons  mitt  onc« 


FIGURE  63.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 
SPECTRUM  L06R2  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION 
Note :  The  plotted  experimental  data  are  in  error  and 
should  be  multiplied  by  a  factor  of  0.50,  resulting  in 
a  peak  transmittance  of  0.0043  at  8.89  |im. 


L23R2  ood  LOWTRAN  (no  o«ro*ol) 


FIGURE  64.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 
SPECTRUM  L23R2  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION 
Note;  The  plotted  experimental  data  are  in  error  and 
should  be  multiplied  by  a  factor  of  3.32,  resulting  in 
a  peak  transmittance  of  0.0199  at  8.96  |xm. 
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transmit  tanc« 


FIGURE  65.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 
SPECTRUM  L23R4  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION. 
Note:  The  plotted  experimental  data  are  in  error  and 
should  be  multiplied  by  a  factor  of  3.28,  resulting  in 
a  peak  transmittance  of  0.0164  at  8.92 /jm. 


L25RJ  and  LOWTRAN  (no  o«ro*ol) 


FIGURE  66.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 
SPECTRUM  L25R3  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION. 
Note :  The  plotted  experimental  data  are  in  error  and 
should  be  multiplied  by  a  factor  of  2.58,  resulting  in 
a  peak  transmittance  of  0.0219  at  8.96  ptm. 
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trontrnlttonc« 


L2SR£  <r»d  LOWTRAN  (no  ooroool) 


FIGURE 


FIGURE 


SPECTRUM  L25R5  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  NO  AEROSOL  ATTENUATION. 
Note:  The  plotted  experimental  data  are  in  error  and 
should  be  multiplied  by  a  factor  of  2.59,  resulting  in 
a  peak  transmittance  of  0.0207  at  8.90pm. 


S25R1  »id  LOWTRAN  (100  km  visibility) 


68.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 

SPECTRUM  S25R1  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  A  VISIBILITY  OF  100  km. 
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S25R4  ond  LOWTRAN  (100  km  visibifty) 


FIGURE  69.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 

SPECTRUM  S25R4  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  A  VISIBILITY  OF  100  km. 


S2bR6  and  LOWTRAN  (K>0  km  viaibility) 


FIGURE  70.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 

SPECTRUM  S25R6  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  A  VISIBILITY  OF  100  km. 
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L2SR5  and  LOttTTRAN  (100  km  vitibUity) 


9 


FIGURE  71.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 
SPECTRUM  L25R5  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  A  VISIBILITY  OF  100  km. 
Note :  The  plotted  experimental  data  are  in  error  and 
should  be  multiplied  by  a  factor  of  2.59,  resulting  iii 
a  peak  transmittance  of  0.0207  at  8.90 /^m. 


7j0  AJ>  9J>  10.0  110  ns>  U.0  U.0 


•ov«l«ngth  (um) 

FIGURE  72.  TRANSMITTANCE  VERSUS  WAVELENGTH  FOR  MEASURED 
SPECTRUM  L25R3  COMPARED  WITH  LOWTRAN  6  CALCU¬ 
LATIONS  PERFORMED  USING  A  VISIBILITY  OF  100  km. 
Note :  The  plotted  experimental  data  are  in  error  and 
should  be  multiplied  by  a  factor  of  2.58,  resulting  in 
a  peak  transmittance  of  0.0219  at  8.96pm. 
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uation  along  the  measurement  path.  The  rationale  behind 
this  omission  is  twofold:  (1)  The  aerosol  attenuation  along 
the  very  long  paths  studied  in  this  experiment  is  expected 
to  be  substantially  smaller  than  the  corresponding  molecular 
absorption  and  the  approach  to  modeling  this  aerosol  atten¬ 
uation  is  quite  uncertain.  Visibilities  estimated  at  sea 
level  during  the  experimental  period  were  nominally  in  the 
range  of  40  km  [17].  However,  outdoor  lights  in  the  vicinity 
of  the  AMOS  observatory  were  clearly  visible  from  the  source 
locations  on  Lanai  over  a  66.2  km  path  using  the  unaided 
eye;  and  (2)  In  the  event  that  transmittance  spectra  cal¬ 
culated  without  an  aerosol  extinction  component  exhibit 
lower  transmittance  values  than  the  measured  spectra,  it  can 
be  clearly  shown  that  the  molecular  absorption  modeling  is 
too  pessimistic. 

This  is  indeed  the  case  for  the  comparisons  shown  in 

Figures  63-67  for  the  7.0-14.0  urn  spectral  band. 

The  comparisons  of  the  measured  transmittance  spectra 
collected  on  25  July  with  the  corresponding  LOWTRAN  calcu¬ 
lations  were  repeated  using  a  value  of  100  km  for  the  visi¬ 
bility  in  the  LOWTRAN  6  Navy  maritime  aerosol  model.  The 
results  of  these  comparisons  are  shown  in  Figures  68-72. 

As  can  be  seen  from  an  examination  of  the  comparisons 
shown  in  Figures  55-62  for  the  3. 0-5.0  urn  band,  the  cal¬ 

culated  peak  transmittances ,  with  no  aerosol  attenuation 
included,  are  typically  about  a  factor  of  2  to  3  greater 

than  the  measured  values.  The  general  shape  of  the  band 
envelope  is  comparable  for  the  measured  and  calculated 

spectra;  however,  the  detailed  structure  of  measured  data  is 
seen  to  vary  for  the  different  measurements.  It  is  felt 
that  the  residual  effects  of  low  frequency  refractive-index 
turbulence,  incompletely  averaged  in  the  measurement,  as 
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well  as  the  intermittent  presence  of  very  thin  clouds  along 
the  line  of  sight  are  responsible,  at  least  part,  for  the 
high  frequency  structure  seen  in  the  experimental  data. 


The  measured  long-wavelength  spectra  shown  in  Figures 
63-67  consistently  show  what  appears  to  be  an  absorption 

feature  occurring  between  9.5  and  10.0  ym.  A  strong  O3 

absorption  band  is  centered  at  9.6  urn  and  might  possibly  be 
responsible  for  this  feature.  No  ozone  absorption  was  in¬ 
cluded  in  the  LOWTRAN  calculation  because  of  the  uncertainty 
in  assigning  a  concentration  for  this  highly  variable  atmo¬ 
spheric  constituent.  One  can  see  from  an  examination  of 

Figures  63-67  that  the  calculated  peak  transmittance  values 
near  9.0  ut\  are  lower  than  the  measured  values  by  about  a 
factor  of  3  to  13.  Again,  no  aerosol  extinction  component  is 
included  in  the  comparisons  shown  in  Figures  63-67. 

In  the  comparisons  shown  in  Figures  68-72  a  value  of 
100  km  was  used  for  the  visibility  in  calculating  the 
aerosol  attenuation  using  the  LOWTRAN  6  maritime  aerosol 

model.  The  general  agreement  in  trend  and  magnitude  for  the 
3-5  vm  comparisons  shown  in  Figures  68-70  is  auite  good. 
The  experimental  transmittances  are  about  10-40%  greater 
than  the  calculated  values,  a  level  of  agreement  which  is 

within  the  range  of  uncertainty  in  the  experimental  data. 


It  is  instructive  to  note  in  these  same  comparisons 
that  the  small  region  of  measureable  transmittance  between 
4.6  and  4.8  urn  is  peaked  at  a  longer  wavelength  in  the 
measured  spectra  than  in  the  calculations.  The  relative 
discrepancies  in  maximum  transmittance  values  in  this  region 
(4. 6-4. 8  wi )  shown  in  Figures  68-70  are  much  larger  (about  a 
factpr  of  2)  than  in  the  region  between  3.3  and  4.3  wn ,  in- 
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dicating  that  the  LOWTRAN  calculation  is  relatively  more 
pessimistic  in  the  long  wavelength  region  of  the  3-5  urn 
atmospheric  window  near  4.6  urn. 

The  7.0-14.0  wi  comparisons  shown  in  Figures  71  and  72 
are  also  based  on  a  value  of  100  km  for  the  visibility  in 
the  LOWTRAN  Navy  maritime  aerosol  model.  Whereas  compar¬ 
ability  in  the  magnitude  of  transmittance  values  is  seen  for 
the  3-5  wm  data  shown  earlier  in  Figures  68-70,.  the  7- 
14  UTi  comparisons  shown  in  Figures  71  and  72  show  that  mea¬ 
sured  transmittance  values  exceed  the  calculations  by  about 
a  factor  of  10. 
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6 

SUMMARY  AND  CONCLUSIONS 


6.1  THE  NASA  DIAL  DATA 

6.1.1  SUMMARY 

Transmi ttances  were  calculated  for  each  of  the  measured 
water  vapor  profiles  for  a  cw  DF  laser  with  the  MIRACL  power 
spectrum,  and  for  a  pulsed  DF  laser  power  spectrum  with  a 
generic  power  spectrum.  All  atmospheric  gases  which  absorb 
in  the  DF  spectral  region  were  included  in  the  calculations. 
With  the  exception  of  water  vapor,  all  gases  were  assumed  to 
be  in  their  midlatitude  summer  concentrations.  The  calcula¬ 
tions  included  line  absorption  and  continuum.  Aerosols  were 
not  included.  Transmlttances  were  also  calculated  for  aver¬ 
age  profiles  based  on  the  data  and  for  standard  atmosphere 
profiles  . 

Thermal  blooming  calculations  were  done  to  determine 
the  importance  of  having  a  shipboard  system  to  monitor  water 
vapor  profiles.  These  calculations  were  also  done  for  both 
cw  and  pulsed  DF  laser  systems. 

These  data  and  the  calculations  based  on  them  represent 
one  day  at  one  particular  location.  They  are  useful  in 
developing  a  perspective  on  the  magnitudes  of  some  of  the 
problems  involved  in  DF  propagation  in  the  real  atmosphere, 
but  sweeping  generalizations  should  not  be  made. 

6.1.2  CONCLUSIONS  AND  RECOMMENDATIONS 

For  the  conditions  represented  by  the  water  vapor  pro¬ 
file  data,  transmittance  for  a  cw  DF  laser  over  a  2-km  path, 
either  vertical  or  slant,  was  about  90%.  For  a  pulsed  DF 
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laser,  the  2-km  transmittance  was  around  85%.  The  pulsed 
transmittance  would  have  been  higher  if  the  long-wavelength 
lines  (4.6  to  4.8  pm)  had  not  been  included.  There  was  only 
a  small  amount  of  fluctuation  about  these  transmittance 
values,  and  they  differed  only  slightly  fron  transmi ttances 
calculated  using  the  standard  atmosphere  water  vapor  pro¬ 
file. 

No  significant  spatial  patterns  were  observed  in  the 
transmi ttances .  In  particular,  no  effect  of  the  gulf  stream 
boundary  on  transmittance  was  observed. 

In  the  thermal  blooming  calculations,  the  performance 
of  a  laser  optimized  using  a  standard  atmosphere  water  vapor 
profile  was  compared  to  the  performance  of  a  laser  optimized 
using  the  profile  generated  from  the  DIAL  data.  Both  lasers 
were  assumed  to  be  operating  in  an  atmosphere  represented  by 
the  DIAL-generated  profile.  The  difference  in  intensity 
delivered  to  the  target  was  small,  less  than  four  percent 
for  both  pulsed  and  cw  systems  and  a  2-km  range.  The  con¬ 
clusion  is  that  on  this  particular  day,  for  this  particular 
scenario,  a  system  to  monitor  the  water  vapor  profile  would 
not  have  been  necessary. 

The  water  vapor  profile  has  only  a  minor  effect  on  DF 
laser  propagation,  either  cw  or  pulsed,  both  in  terms  of 
linear  and  of  nonlinear  degradation  of  the  beam.  Based  on 
these  data,  precise  knowledge  of  the  water  vapor  profile 
does  not  appear  to  be  critical. 

Detailed  treatment  of  the  propagation  of  DF  laser  ra¬ 
diation  can  be  quite  complicated,  as  a  high-energy  beam  in 
this  spectral  region  can  have  some  complex  nonlinear  inter¬ 
actions  with  atmospheric  N2 ,  ^2©,  CO2 ,  and  aerosols.  Of 

this  group,  only  aerosols  are  highly  variable.  Aerosols 
were  not  considered  at  all  in  this  study.  Aerosol  data  were 
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not  supplied,  and  aerosol  models  were  not  included  in  the 
transmittance  modeling  or  the  thermal  blooming  calculations. 
Future  work  on  the  effects  of  atmospheric  variability  on  DF 
laser  radiation  should  probably  concentrate  on  aerosols 
rather  than  water  vapor.  Open  ocean  data  are  needed  to 
verify  existing  aerosol  models  or  for  use  in  generation  of 
better  models.  Nonlinear  effects,  especially  aerosol  en¬ 
hancement  of  thermal  blooming,  may  be  important  in  DF  propa¬ 
gation,  and  study  of  this  subject  reguires  reliable  informa¬ 
tion  about  the  aerosols  encountered  at  sea.  The  importance 
of  aerosol  growth  with  humidity,  and  the  resulting  depen¬ 
dence  of  the  aerosol  profile  on  the  water  vapor  profile, 
should  be  established;  the  DIAL  system  seems  uniquely  suited 
for  measurements  related  to  this  subject. 


6.2  THE  LFLID  SUPPORT  DATA 

The  long-path  atmospheric  transmittance  and  radiosonde 
measurements  performed  in  support  of  the  LFLID  tests  were 
carried  out  during  the  month  of  July  1984.  After  editing 
the  transmittance  data  to  remove  those  measurements  which 
were  affected  by  cloud  occurrences  along  the  measurement 
line  of  sight,  data  collected  during  the  early  morning  hours 
( HST )  on  four  days  were  selected  for  further  analysis.  The 
measurement  dates  and  times  for  these  data  are  given  in 
Table  5,  Section  2.2.  Figures  15  and  16  in  Section  2.2  show 
the  vertical  temperature  and  water  vapor  profiles  determined 
from  radiosonde  data  which  were  collected  during  the  four 
selected  measurement  periods. 

The  temperature  profiles  show  the  presence  of  inver¬ 
sions  occurring  at  altitudes  between  1.5  and  2.1  km  during 
three  of  the  four  periods.  Reduced  water  vapor  densities 
above  1.8  km  appear  to  be  correlated  with  the  temperature 
inversions . 
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The  results  of  the  long-path  transmittance  measurements 
and  corresponding  LOWTRAN  calculations  are  shown  in  Figures 
55-72  in  Chapter  5.  An  examination  of  the  measured  and  cal¬ 
culated  3-5  wn  transmittance  comparisons  shown  in  these 
figures  indicates  that  reasonably  good  agreement  between  the 
measured  data  and  the  LOWTRAN  6  model  exists.  Two  variations 
of  the  LOWTRAN  calculations  were  performed,  one  omitting  any 
aerosol  attenuation  contribution  whatsoever,  and  the  other 
using  a  nominal  ‘visibility  of  100  km.  It  can  be  seen  from 
an  examination  of  the  plotted  results  shown  in  Figures  60- 
62,  compared  with  those  shown  in  Figures  68-70,  that  the  two 
LOWTRAN  calculations  bracket  the  experimental  data  in  the  3- 
5  un  region.  The  100  km  visibility  calculation  is  seen  to 
agree  with  the  measured  data  to  within  the  estimated  range 
of  experimental  uncertainty. 


Comparisons  similar  to  those  in  the  3-5  wn  region  were 
carried  out  for  the  long-wavelength  (7.0-14.0  ym)  band  and 
are  shown  in  Figures  63-67  and  Figures  71  and  72.*  The  cal¬ 
culations  performed  using  a  visibility  of  100  km  show  that 
the  LOWTRAN  model  underestimates  the  measured  transmittance 
values  by  about  a  factor  of  ten  (Figures  71  and  72)  while 
the  comparable  comparisons  in  the  3-5  urn  region  (Figures  68- 
70)  show  good  agreement. 
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This  observation  is  consistent  with  results  obtained  in 
an  analysis  of  atmospheric  CO2  laser  transmittance  mea¬ 
surements  in  the  9-11  wn  region  performed  by  NRL  in  previous 
experiments  at  Cape  Canaveral,  Florida  [18]. 

The  experimental  data  presented  in  Figures  68-70  show 
that  the  measured  transmittance  in  the  region  between  4.6 
and  4.9  un  is  about  a  factor  of  two  larger  than  the  calcu¬ 
lated  value,  indicating  that  the  model  is  relatively  pessi¬ 
mistic  in  this  sub-region  when  compared  with  the  good  agree¬ 
ment  evident  in  the  3. 3-4. 3  wn  sub-region. 

The  measured  transmittances  in  the  7-14  wm  region  con¬ 
sistently  show  the  presence  of  an  absorption  feature  between 
9.5  and  10.2  rni.  This  feature  is  very  likely  due  to  the 
presence  of  ozone  in  the  propagation  path.  The  uncertainty 
associated  with  determination  of  an  appropriate  O3  concen¬ 
tration  for  use  in  the  LOWTRAN  calculations  precluded  incor¬ 
poration  of  the  O3  absorption  features  into  the  calculated 
transmittance  spectra. 
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